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MELTED  FLUX  LIQUIDS  IN  HIGH -Tc  SUPERCONDUCTORS 


David  R.  NELSON 

Lyman  Laboratory  of  Physics,  Harvard  University,  Cambridge,  Massachusetts  02138 


A  theory  of  the  entangled  flux  liquids  which  arise  in  the  new  high-Te  superconductors  is  reviewed.  New  physic* 
appears  because  of  the  weak  interplanar  couplings  and  high  critical  temperatures  in  these  materials.  Flux  line 
wandering  melts  the  conventional  Abrikosov  flux  lattice  over  large  portions  of  the  phase  diagram  and  leads  to 
a  novel  entangled  vortex  state.  We  suggest  that  a  heavily  entangled  flux  liquid  could  exhibit  glassy  behavior 
on  experimental  time  scales,  in  analogy  with  viscoelastic  behavior  in  dense  polymer  melt*. 


1.  INTRODUCTION 

One  of  the  many  fascinations  of  the  CuO}-based 
superconductors^1)  is  the  possibility  of  novel  fluctua¬ 
tion  effects  due  to  the  high  critical  temperatures  and 
small  coherence  lengths.  Although  fluctuations  are  usu¬ 
ally  limited  to  the  immediate  vicinity  of  critical  point*, 
it  was  recently  argued  that  fluctuation*  in  the  high-Te 
materials  lead  to  a  new  entangled  flux  liquid  pKatt  in 
a  magnetic  field,  due  to  flux  line  wandering  as  vortex 
filament*  traverse  the  sample.1’* 

These  speculations  are  consistent  with  striking  vi¬ 
brating  reed  experiments  by  Gammel  et  of.,4  who  find 
signal  suggestive  of  flux  lattice  melting  in  single  crys¬ 
tal*  of  B^SrjCaCujOi  (BSCCO)  at  temperatures  well 
below  the  onset  of  the  Meissner  effect. 

Melted  flux  liquids  are  already  familiar  from  discus¬ 
sions  of  two-dimensional  superconducting  films:  Dislo¬ 
cation-mediated  melting  of  the  flux  lattice  leading  to 
both  ordinary  and  hexatic  liquid  phase*  of  essentially 
point  vortices  was  explored  theoretically  several  years 
ago  by  D.  S.  Fi*her.*  The  novelty  of  high-T,  super¬ 
conductor*  lie*  in  the  possibility  of  a  melted  liquid  of 
entangled  tin*  defect*  in  three  dimension*. 

We  can  estimate  when  the  Abrikosov  theory  break* 
down  from  a  simple  random  walk  argument.1  We  con¬ 
sider  a  single  flux  line  f[i)  and  determine  how  far  it 
wanders  perpendicular  to  the  *  axis  (fi  ||  i)  as  it  tra¬ 
verses  the  sample.  In  Ref.  3,  it  is  shown  that  the  vortex 
“diffuse**  a*  a  function  of  the  time-like  variable  *, 


v/<irt*)  -rlO)!1)  =  y/wi  , 

where  the  “diffusion  constant*  D  is  given  in  term*  of 
the  mass  anisotropy  M^/M,  a  10-1  and  the  vortex 
line  tension  ej.by 

n  -  b. 'I 

M.  e;  ' 

At  T  =  77  K,  we  take  Hcl  a  101  G  (which  deter¬ 
mines  e j )  and  <s  101  and  find  D  =  10~*  cm, 

so  that  vortex  lines  wander  a  distance  of  order  1  micron 
while  traversing  a  sample  of  thickness  0.01  cm.  This 
diffusive  behavior  implies  an  entanglement  correlation 
length,1,1  (*  =  1/2 Dn,  which  is  the  distance  required 
to  wander  a  vortex  spacing  in  a  vortex  liquid  with  areal 
density  n  =  B/6 o-  Collisions  and  entanglement  of  vor¬ 
tex  lines  will  significantly  alter  the  Abrikosov  theory 
whenever 

L»(,  , 

where  L  is  the  sample  thickness. 

2.  LINDEMANN  CRITERION  AND  PHASE 
DIAGRAM 

To  determine  where  the  flux  lattice  melts,  we  as¬ 
sume  that  the  externa]  field  is  aligned  with  the  z- 
direction,  perpendicular  to  the  CuOj  planes,  and  de¬ 
scribe  the  trajectory  of  the  i-th  vortex  by  a  function 
r,(z).  If  the  average  position  of  the  i-th  vortex  in  a 


triangular  lattice  denoted  by  £,,  we  can  define  a  two- 
dimensional  displacement  field  u(/Z,,z)  by 


r~.(2 


-  «(£,  a) 


In  the  continuum  limit,  this  displacement  becomes  a 
function  u  =  u(x,y,  z).  The  excess  free  energy  *G{u(r)} 
associated  with  small  gradients  of  a  is4 


\  j  d*r  +  U]s  -t-  K  (^j 


where 


,  -  1  /  du„  dua\ 

0  =  H^7  + ’  a’/3  =  I'v 


is  the  symmetrised  two-dimensional  strain  matrix,  ^ 
and  A  are  Lam*  coefficients,  and  K  is  a  till  elastic 
constant.  The  mean-square  displacement  of  u  vortex  is 
then  (assuming  2n  +  A  >>  *») 


(l“(r'o)l*) 


[EhL 

V  4* 


According  to  the  Lindemann  criterion,  melting  wifi  oc¬ 
cur  when  y/(u*)  exceeds,  say,  ^  of  the  lattice  con¬ 
stant.  To  evaluate  Eq.  (3),  we  need  estimates  of  n  and 
K. 

Although  the  shear  modulus  will  be  of  order  its 
value  in  isotropic  superconductors/  the  tilt  modulus 
K  will  typically  be  much  tmaller  than  the  isotropic 
value®  K  =  when  the  anisotropy  of  high-Ti  ma¬ 
terials  is  taken  into  account.  A  reasonable  guess,  valid 
for  ffct  «  B  <<  Bn,  is* 


Afj.  BB 
Afi  tr 


(1) 


where  Mj.  is  the  in-plane  effective  mass,  and  M,  a 
105Afj_  is  the  much  larger  mass  describing  the  weak 
coupling  between  CuOj  planes.  The  tilt  modulus  is 
smsll  because  the  magnetic  field  does  not  rotate  appre 


ciably  when  vortex  lines  are  tilted  in  a  highly  anisotropic 
superconductor.*  The  field  remains  essentially  uniform 
and  perpendicular  to  the  CuOj  planes  whenever  d  « 
A,  where  d  is  the  vortex  spacing  and  A  is  the  London 
penetration  depth.  The  estimate  (4)  breaks  down  near 
Bo  because  most  of  the  material  is  in  its  normal  state, 
and  the  macroscopic  magnetic  properties  must  then  re¬ 
flect  the  (small)  magnetic  anisotropy  of  a  normal  metal. 
As  pointed  out  by  D.  Huse,®  Eq.  (1)  also  breaks  down 
near  Bei,  because  of  the  magnetic  dipole  coupling  be¬ 
tween  vortices  present  for  d  >  A.  Over  the  large  range 
of  fields  where  (i)  is  valid,  it  is  easy  to  estimate  the  root 
mean  square  phonon  displacements  and  show  that  the 
BSSCO  fiux  lattice  is  indeed  melted  over  a  wide  range 
of  fields  and  temperatures.* 

The  Lindemann  criterion  assumes  wave  vector-inde- 
per  J*nt  elastic  constants  which  can  be  a  gros*  oversim¬ 
plification,  especially  if  fluctuations  are  important.  A 
more  accurate  theory  would  use  wave  vector  dependent 
elastic  constants.*  This  has  been  done  by  Houghton, 
Pelcovits,  and  Sudbo,10  who  generalize  nonlocal  elas¬ 
tic  constants  derived  by  Brandt.7  As  pointed  out  by 
Houghton  et  al.,  there  are  significant  corrections  to  the 
Lindemaan  criterion  due  to  nonlocal  effects  in  high-Te 
materials.  These  nonlocal  effects  act  to  enhance  the 
fluctuation  induced  melting  considt>ed  above.  By  using 
the  Lindemann  ratio  as  a  fitting  parameter  (they  find 
vV)  =  0.4d  at  melting),  these  authors  are  able  to 
obtain  good  fits  to  the  melting  curves  of  Ref.  6  for  both 
YBCO  and  BSCCO. 

The  uWemann  analysis  sketched  above  cannot  be 
correct  in  the  vicinity  of  Bcl,  i.e.,  when  the  flux  line 
spacing  is  of  order  the  London  penetration  depth  or 
more.  This  is  the  regime  of  most  flux  decoration  experi¬ 
ments.  The  classical  theory,  based  on  energetic  consid¬ 
erations  alone,  predicts  that  n  ~~  A  ~  exp(-d/A)  as  the 
vortex  separation  d  -»  oo.7  There  are,  however,  im¬ 
portant  entropic  effects  due  to  flux  line  wandering  which 
dominate  over  flux  line  interaction  energies  in  the  dilute 


limit.  It  can  be  shown  that  the  effect  of  these  fluctu¬ 
ation*  is  to  produce  a  thin  sliver  of  melted  flu*  liquid 
phase  ju*t  above 

Figure  1  shows  a  schematic  phase  diagram  for 
BSCCO  which  combines  the  high  field  Lindemann  cri¬ 
terion  melting  curves  of  Houghton  e<  of.10  with  the  low 
field  estimate  of  the  melting  curve  discussed 


Fig.  1.  Schematic  phase  diagram  for  BSSCO  with  H 
directed  perpendicular  to  the  CuOj  planes. 

above.11,13  The  field  is  assumed  to  be  perpendicular  to 
the  CuOj  planes.  The  region  occupied  by  the  entan¬ 
gled  flux  liquid  is  smaller  for  YBCO  than  for  BSCCO, 
because  M,/M±  is  much  smaller  in  this  case.  For 
most  conventional  superconductors,  the  reentrant  melt¬ 
ing  curve  H,(T)  is  indistinguishabiy  close  to  He\(T) 
and  Hci(T).  The  mean-field  transition  line  at  £fcl 
need  not  be  a  sharp  phase  transition  when  fluctuations 
are  taken  into  account; ,3’1*  experimentally  it  is  defined 
by  the  onset  of  an  (incomplete)  Meissner  effect  with  de¬ 
creasing  temperatures. 

3.  PROPERTIES  OF  THE  ENTANGLED  FLUX 
The  entangled  flux  liquid  is  a  new  phase  of  matter. 


stabilised  by  the  entropy  available  when  flux  lines  en¬ 
tangle.  Of  particular  interest  is  the  vot  ex  line  corre¬ 
lation  function,  which  can  be  measured  with  neutron 
scattering.14  These  experiment*  measure  fluctuations  in 
the  local  density  of  flux  lines,  defined  by 

N 

=  YL  ~  • 

i=i 

Neutron  scattering  measures 

where  n(q :±,qM)  is  the  Fourier  transform  of  n(r"x,r). 
The  physics,  however,  is  more  conveniently  discussed  in 
terms  of  the  partial  Fourier  transform 

S(9i,z)=  r 

j -oo  ( 

=  (n(?i,r  +  eo)nl{9j.,zo))  ■  (2) 

This  con-elation  function  reduces  to  the  structure 
function  of  a  cross-section  of  the  vortex  lines  when  x  =  j 
0.  The  correlations  in  any  such  con- '.ant  x  cross- 
section  should  be  similar  to  those  of  a  two-dimensional 
liquid.  Equation  (2)  describes  more  generally  how  the  i 
Fourier  components  of  this  2d  structure  function  de¬ 
cay  along  the  z-axis  due  to  entanglement.  As  shown 
in  Ref.  3  (by  exploiting  an  analogy  with  superfluidity 
in  two-dimensional  bosons)  this  decay  can  be  approxi¬ 
mated  by 

5(#l,e)»5(,i,s  =  0)e-(^)''/‘»r  ,  (3) 

where  e(gx)  is  the  phonon-roton  excitation  spectrum  of 
the  corresponding  superfiuid.  The  excitation  spectrum 
defines  via  Eq.  (3)  a  ^-depend  ent  correlation  length 

£||(?i)  =  *BT/e(qL)  . 

The  quantity  f;;(?i)  equals  the  ‘entanglement  corre¬ 
lation  length*  (,  when  ~  1/d,  where  d  is  an 


intervortex  spacing.* 

There  is  also  interesting  modifications  of  the  B{R) 
constitutive  relation  near  Etei  where  the  transition  from  i 
the  Meissner  phase  is  into  an  entangled  flux  liquid. 3,3 
A  discussion  of  the  hydrodynamics  of  equilibrated  flux 
liquids  is  in  preparation. IS 

4.  SLUGGISH  DYNAMICS  DUE  TO  ENTANGLE- 
MENT 

Disorder  will  be  most  effective  in  suppressing  flux  , 
flow  resistivity  in  vortex  phases  which  exhibit  a  nonzero 
shear  modulus  on  experimental  time  scales.  In  con¬ 
ventional  superconductors,  for  example,  a  few  securely 
pinned  flux  lines  can  pin  the  entire  flux  lattice  if  there  is 
large  shear  rigidity.  We  argue  here  that  very  long  shear 
relaxation  times  may  also  be  possible  in  an  entangled 
flux  liquid.*  Although  it  will  greatly  affect  the  response 
to  a  dilute  concentration  of  pinning  centers,  this  mech¬ 
anism  for  obtaining  a  shear  modulus  does  not  require 
pinning  centers  to  work. 

Figure  2  shows  a  heavily  entangled  flux  liquid,  view-  ; 


dimensional  random  walk,  occasionally  entwining  around 
one  if  its  near  neighbors.  Figure  2  is  like  a  two-dimension¬ 
al  polymer  melt,  with  the  understand!)  g  that  polymer* 
which  appear  to  cross  in  ;ne  projection  are  actually  or¬ 
dered,  one  above  the  ather,  according  to  which  of  the 
crossing  lines  has  the  largest  *  coordinate.  The  s 
coordinate  plays  the  role  of  a  monomer  index  in  this 
polymer  analogy. 

If  the  configur»tion  of  vortex  lines  in  Figure  2  is 
subjected  to  a  shear  stress  oty,  the  constraints  pro¬ 
vided  by  the  entanglements  in  this  figure  will  make  re¬ 
laxation  difficult.  The  notion  of  constraints  requires  a 
large  barrier  to  flux  line  cutting.  When  two  lines,  both 
approximately  parallel  to  3,  cross,  they  must  pass  over 
a  barrier  of  order  2<i/,  where  ti  is  the  line  tension, 
and  1  is  the  distance  along  f}  over  which  the  crossing 
occur*.  This  barrier  arises  because,  as  two  lines  cross, 
they  behave  like  single  double  quantized  vortex  over  the 
distance  l.  The  supercurrents  responsible  for  this  bar¬ 
rier  flow  primarily  within  the  CuOj  planes,  due  to  the 
anisotropy.*  Upon  expressing  Cj  in  terms  of  the  lower 
critical  field,  we  find 


2tif 

kBT 


Hc\M 

2xkgT 


ss  50  , 


(<) 


where  we  ha«e  set  Hc i  =  80  Oe,  l  =  10A,  and  T  - 
77  K.  Thus,  line  crossings  will  be  very  difficult,  just  as 
for  polymer*  in  an  entangled  melt. 

The  analogy  with  polymer  melts'*  suggest  that, 
after  an  initial  transient,  the  strain  uIy  will  be  pro¬ 
portional  to  the  applied  stress  for  times  t  <  rL, 


U x y  SS  <T»y/p O  i 

which  defines  a  shear  modulus  for  times  shorter 
than  the  relaxation  time  rj,.  For  times  longer  than  rj., 
the  material  will  behave  like  a  viscous  liquid,  with  the 
strain  rising  linearly  with  (, 


Fig.  2.  Schematic  top  view  of  an  entangled  flux  liquid, 
ed  down  the  s  axis.  Each  flux  line  executes  a  two- 


ttx*  a*  Oxyth  • 


which  define*  *  shear  viscosity  f).  The  key  quantity 
in  this  description  of  the  viscoelastic  response  is  the 
relaxation  time  rif  which,  for  polymers,  is  a  strong 
function  of  the  polymerisation  index. 

The  analogous  time  scale  in  entangled  flux  liquids 
will  be  a  sensitive  function  of  the  sample  thickness  L. 
To  estimate  this  time,  we  use  the  de  Gennes  reptation  ' 
theory  of  polymeric  relaxation  in  three  dimension*,1®  in 
which  each  polymer  is  assumed  to  diffuse  along  a  “tube* 
defined  by  the  entanglement  constraints.  A  very  similar 
situation  arises  in  the  two-dimensional  polymeric  prob-  | 
lem  considered  here.  The  de  Gennes  theory  predicts 
that  ti  varies  approximately  as  the  cube  of  the  sam-  j 
pie  thickness 

a  =  ML/a* . 

where  tq  is  a  microscopic  time,  and  we  expect  that  £' 
is  of  the  order  of  the  entanglement  correlation  length. 
Motion  of  tine  center  of  mass  of  a  single  vortex  line  is 
described  by  the  diffusion  constant18 

Dl  —  D0(L/i'a)~i  , 

where  D0  is  the  diffusion  constant  of  a  point  vortex 
in.  say,  an  isolated  CuOj  plane.  A  rough  estimate  is 
D0  ~  fc/m,  —  1  cmJ/sec.  The  microscopic  relaxation 
time  should  be  of  order 

Tq  ~  \  J Don  ~  ^0/ DqB 

If  the  intervortex  spacing  is  d  —  {60/ B)1?1  a  100  A  and 
{|  =  2000  A,  we  find  from  Eq.  (13)  that  =  103  sec 
or  about  15  min  for  L  =  1  cm.  The  shear  viscosity  in 
this  regime  should  be  of  order18 

n  «  riMO  , 

or  i)»2x  10*  poise  using  standard  estimates3  of  po 
with  H  =  \Hcj. 

The  ideas  sketched  above  suggest  that  regimes  of 
very  high  viscosity  (over  10*  times  that  of  water)  are 


possible  in  sufficiently  thick  samples.  Because  the  line 
tension  (1  vanishes  near  Hc i,  one  must  be  sufficiently 
far  below  the  mean  field  Hct(T)  before  the  barrier  (4) 
associated  with  entanglement  becomes  effective.  One 
might  even  speculate  that  the  sluggish  dynamics  which 
results  with  decreasing  temperature  is  the  origin  of  the 
“irreversibility  lines*  like  those  discussed  by  Malotemoff 
it  aJ.‘7 
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ANISOTROPY  IN  MAGNETIC  PROPERTIES  OF  HIGH-Tc  MATERIALS 
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A  brief  review  is  given  on  recent  results  concerning  equilibrium  flux-line  lattices,  the  field 
distribution  within  a  cell  of  vortex  lattices,  the  transverse  magnetization  and  the  anisotropy 
related  torque.  Current  carrying  properties  of  boundaries  between  misaligned  anisotropic 


superconductors  are  discussed. 

This  brief  review  is  intended  to 
summarize  the  recent  results  and  to  provide 
a  reference,  however  selective  and 
incomplete,  for  anyone  entering  the  rapidly 
expanding  subject  of  the  title.  The  simple 
way  to  describe  major  anisotropy  effects  is 
to  replace  in  the  London  or  Ginzburg-Landau 
equations  the  squares  of  the  isotropic 
penetration  depth  and  of  the  coherence 
length,  X2  and  l;2,  with  tensors  (X2)ik  = 
/i-mik  and  (^)jk  =  ^m^-l,  respectively. 
Tensor  mjk  (usually  called  "mass-tensor", 
although  it  is  related  to  the  superfluid 
density  and  incorporates  all  sources  of 
anisotropy  of  the  tensor  type)  is  normalized 
on  the  average  mass  (M]M2M3)1/3,  where 
M's  are  the  principal  values  of  Mjk.  The 
average  values,  X  and  are  related  to  the 
principal  values:  X^XjXjXy,  Note, 

that  only  the  "squares",  (X2);^  and 
transform  as  tensors  under  coordinate 
transformations.  Also,  not  all  anisotropy 
effects  can  be  described  with  this  simru 
approach;  e  g.,  a  possible  anisotropy  in  the 
basal  plane  of  tetragonal  materials  is 
completely  missed  in  this  description.  Still, 
in  materials  of  interest  the  basal  plane 
anisotropy  is  usually  weak,  and  m;^ 

approach  provides  a  simple  and  unified 
description  of  various  effects  of  the  strong 


"c-a"  anisotropy. 

The  only  restriction  for  the  London 
theory  to  hold  is  a  spaciaily  constant 
moduli)'  of  the  order  parameter.  This  is 
fulfilled  in  Fields  H«HC2.  where  the  vortex 

cores  of  the  size  4  are  not  overlapped,  the 
usual  Situation  for  high-Tc  materials.  The 

absence  of  any  restriction  upon  T  and  the 
simplicity  make  the  London  approach  a 
useful  tool  in  studying  magnetic  properties 
of  high-Tc's.  We  will  concentrate  in  the 

following  on  results  obtained  within  this 
approach. 

Using  the  London  free  energy1 -2 

,  2  2 

F=J(h  +  X  mucurl^curlkh)  dv/8)t  ,  ( 1 ) 

one  can  compare  different  flux  arrange¬ 
ments  and  find  one  corresponding  to  the 
thermodynamic  equilibrium.  For  inter¬ 
mediate  applied  Fields  Ha,  «  Ha«Hc2, 
this  was  done  in  Ref.  3  for  the  uniaxial  case. 
The  magnetization  M  in  this  domain  is  on  the 
order  of  Hcj  and  therefore  is  small  with 
respect  to  B,  H  and  Ha;  the  three  latter 
vectors  are  almost  the  same  :r  '"*h>e  and 
direction.  This  simpliFies  the  problem 
allowi;  g  to  disregard  demagnetization 
effects.  For  the  Field  at  an  angle  8  to  the  c 
direction  of  the  crystal  the  flux  line  lattice 
(FLL)  in  the  plane  normal  to  B  is  made  of 
isosceles  triangles  with  the  base-to-side 


ratio  of  V  l+3mym(8)  12  where  m(8)  = 
masin-8  +  mccos- j  The  short  base  is  directed 
along  the  r  -lion  of  c  onto  the  plane 
normal  'u  b.  the  long  height  of  the  triangle  is 
sainted  in  the  ab  plane.  The  FLL  t s  locked 
upon  the  crystal ;  any  rotation  of  the  FLL  as  a 
whole  round  B  causes  an  energy  increase, 
the  feature  absent  in  isotropic  materials. 

The  same  equilibrium  FLL  has  been  found 
near  Hcj,2  near  Hc-i.4  and  as  a  particular 
case  in  a  more  general  context  of  3  different 
masses  5  Decoration  experiments  for  Bwc  on 
Y-123  are  consistent  with  London 
predictions. ^ 

Given  the  equilibrium  FLL.  one  can 
study  the  Meld  ,.(ri  within  an  FLL  cell.  This 
is  done  by  soNing  the  London  eqs  in  the 
Fourier  space  and  then  by  the  numerical 
reconstruction  of  htr).  Three  relevant 
results  are  ,o  be  mentioned: 

(aj  For  an  arbitrary  orientation  of 
vortex  axes  (of  B;  within  a  crystal,  the  local 
hir)  has  a  non-zero  component,  h,r(r), 

transverse  to  B.  Thus,  h  is  not  parallel  to 
B.  unless  B  points  in  one  of  the  principal 
directions.  Physically,  hjj.  is  related  to  the 

tendency  of  persistant  currents  to  flow  along 
directions  with  small  masses  which  would 
reduce  the  kinetic  contribution  to  the  free 
energy  (the  second  term  in  Eq.(l)). 

(b)  The  transverse  field  compars  in 
value  to  the  familiar  longitudinal  field 
h  i  g  =  h  B/B,  provided  the  average  vortex 
separation  exceeds  X  (near  Hcj).'  Infields 
well  above  Hc[,  hg.  is  on  the  order  Hcj,  i.e. 

it  is  small  with  respect  to  However,  in 

any  field,  htr  scales  with  the  variable  part 

of  h\g.  Roughly,  htr  scales  with  the 
macroscopic  magnetization  M;  htr  cannot  be 

neglected  as  long  as  properties  related  to  the 
field  variation  (such  as  M)  are  concerned. 


(c)  The  average  <htr>  over  the  FLL  cell 
is  zero.  This  manifestation  of  the  flux 
quantization  has  been  established  in  Ref.  9 
for  H  close  to  Hc2-  The  argument  used  holds 

in  any  field.  Numerical  results  of  Ref.  8 
comply  well  with  this  restriction. 

Due  to  the  "extra"  energy  dependence 
upon  rotations,  the  elastic  properties  of  FLL  s 
in  anisotropic  superconductors  are  richer 
than  in  the  isotropic  case.  Second  order 
terms  of  the  energy  expansion  in  derivatives 
of  the  displacement  u.  contain  not  only 
terms  quadratic  in  strains  u,^  = 

(du^dxk  +  du^/dxj) /2,  but  also  those  quadratic 
in  rotations  ojjj.  =  (duj/dxy.  -  bu^ibx^t!2,  as 
well  as  ones  proportional  to  Uj^toju.  The  2D 
elastic  constants  for  the  FLL  deformations 
that  conserve  the  direction  of  vortex  axes, 
has  been  considered  in  Ref.  10.  A 
potentially  relevant  conclusion:  the 
anisotropy  influence  upon  the  elastic  pro¬ 
perties  of  FLL's  is  much  stronger  than  upon 
such  common  features  of  interest  as  Hc2- 
Hc[  or  M.  It  turned  out,  e.g„  that  in  2D 
there  exist  two  independent  shear  modulae, 
^-66. hard  ^66, easy  w^ose  ratio  is 
(mc/m(0))2  (which  amounts  to  (mc/ma)^  for 
Blc).  Recall  that  the  ratio  of  Hc2's  in  two 
principal  directions  is  only  (mc/ma)^“.  The 
hard  shear  mode  corresponds  to 
displacements  along  the  base  of  the 
elementary  triangle,  whereas  the  easy  one 
describes  the  shear  defor-mation,  in  which 
displacements  are  in  the  ab  plane  (along  the 
long  dimension  of  the  trianglel.  For  the  field 
in  the  ab  plane  even  the  "moderately" 
anisotropic  Y-123,  (mc/ma)^“=5,  the. ratio 
of  the  two  shear  modulae  exceeds  600.  In 
the  Bi-2212  the  estimates  of  the 
(mc/ma)l/2  range  from  17  t Ref.  16)  to  more 
than  50  (Ref.  19'.  putting  the  ratio  of  shears 


between  105  and  102.  The  FLL,  therefore,  is 
very  stiff  with  espect  to  the  hard  shear 
mod'*  and  is  very  soft  with  respect  to  the 
easy  one  tin  which  vortices  slide  along  the 
ab  plane).  The  decoration  results  of  Dolan  ei 
at  ^  for  B^e  are  consistent  with  the  lattice 
being  very  soft  with  respect  to  "sliding" 
deformations.  On  the  other  hand.  Biggs  et 
at  concluded  that  there  is  no  observable 
creep  of  vortices  (Bic)  along  c  in  8i-22I2 
crystals.  This  is  expected  if  the  corres¬ 
ponding  shear  modulus  (Cb6 .hard  'n  ^'s 
case)  is  large.  The  failure  of  the  FLL  (again 
for  Blc)  to  show  a  detectable  response  to 
the  Lorentz  force  along  c  in  experiments 
with  oriented  films  of  Bi-2212  (Ref. 12)  and 
of  Tl-2212  (Ref. 13)  may  also  be  interpreted 
as  an  indication  of  a  large  C^.bard  (,n  *  ^ ' 
123  single  crystal  this  response  is  found  to 
be  finite  but  weak.  Ref.14). 

One  of  the  macrosopic  effects  of  the 
non-zero  transverse  field  in  vortices  of 
anisotropic  materials  is  the  transverse  (to  B) 
component  of  the  magnetization  M  tr.  One 

can  evaluate  \1  if  the  free  energy  of  the  FLL 
i  is  known  as  a  function  of  B.^-l^  The  M[r 

can  be  measured  directly' ^  or  by  observing 
the  torque  T=  IMxBI  =  MtrB.  In  intermediate 

fields  the  torque  data'2  for  Y-123  are  in  a 

1  O 

good  agreement  with  the  London  theory.10 
The  data  should  be  taken  at  sufficiently  high 
temperatures  and  fields,  where  M  is  revers¬ 
ible  to  assure  that  the  torque  measured  is  an 
intrinsic  anisotropy-related  property  (as 
opposed  to  torques  due  to  a  trapped  flux  or 
due  to  the  shape  effects).  The  low-Jc 

material  Bi-2212  provides  a  broad  window 
in  the  T-H  plane  where  such  measurements 
can  be  done;  the  ratio  (n^/mj)1^  extracted 
from  the  torque  data  exceeds  50. *  ^ 

It  should  be  pointed  out  that  among 
various  methods  of  extracting  the  anisotropy 


parameter  y=vmpm,  ,  the  Hc]  measurements 
show  consistently  lower  values  than  those 
obtained  by  other  means.  For  Y-123 
estimates  of  Hcj  produce  y-3,20  while 
magnetization. torque,'  '  Hc2.  and 
decoration^  yield  y-5  or  more.  The 
discrepancy  is  pronounced  in  the  case  of  Bi- 
2212.  for  which  the  Hc  j  data  suggest  again 

v-3  22.23  w;,creas  the  torque  technique 
generates  y-55, '  ^  and  the  H~2  gives  about 
6C  24  The  sut,ject  of  the  Hcj  anisotropy 

deserves  further  study.  At  this  stage,  one 
may  speculate  that  the  maximum  Hcj  for 

H  lie  might  be  underestimated  if  the  crystal 
has  planar  defects  parallel  to  c.  which  allow 
the  vortices  to  penetrate  the  sample  at 
H<HC].  The  minimum  Hc;  for  Hlc.  on  the 

other  hand,  can  be  over-estimated  due  to 
the  platelet-like  crystal  shapes  and  possible 
surface  barriers  in  fields  parallel  to  platelet 
faces.  A  small  misalignment  between  the 
applied  field  and  the  ab-plane  may  again 
result  in  an  overestimate  of  the  ”HC]  parallel 

to  ab".  While  it  is  conceivable  that  the  first 
of  the  reasons  given,  is  responsible  for  an 
under-estimate  of  Hcj  max/Hcl  .min  for  Y- 
123,  it  is  not  yet  clear  whether  the  last  two 
are  sufficient  for  resolving  the  discrepancy 
for  the  Bi-compotnd. 

In  order  to  address  problems  of  the 
current  flow  in  polycrystalline  samples 
within  the  London  approach,  one  has  to 
formulate  the  conditions  at  the  boundaries 
between  misaligned  anisotropic  super¬ 
conductors.25  Perhaps,  the  simplest  way  to 
obtain  these  conditions  has  been  suggested 
by  N.Schopohl:  integrate  the  London  eqs., 
h-t-X2curlI=0,  where  the  vector  Ijj  = 
(4tt/c)mknjn,  over  the  area  of  a  contour 
which  crosses  the  interface.  When  the 
contour  shrinks  in  a  direction  normal  to  the 


boundary  so  that  its  long  sides  run  along  the 
interface  in  the  opposite  directions,  the 
integral  of  h  vanishes.  Thus, 
j  "url  I-ds  =  I-dl  =  dt  i.e.  the  tangential 

components  of  I  must  be  the  same  on  both 
sides  of  the  interface:  t mtnj n ) ieft  = 

*  mtn-in)right  •  t  denotes  any  of  the 

tangential  directions.  This  peculiar  condi¬ 
tion  means  that  for  a  general  misalignment 
where  on  both  sides  may  have  non-zero 

off-diagonal  components,  the  langen-tiai  and 
normal  components  of  j  at  the  interface  are 
coupled.  In  general,  any  transport  current 
through  the  boundary  must  be  accompanied 
by  a  _ urrent  along  the  boundary.  The  latter 
flows  in  opposite  directions  at  opposite  sides 
of  the  interface  and  creates  an  "extra" 
magnetic  field  parallel  to  the  interface,  that 
decays  with  the  distance  from  the  boundary. 
As  a  result,  boundaries  between  anisotropic 
superconductors  may  serve  as  channels, 
along  which  vortices  can  penetrate  the 
sample  at  a  field  smaller  than  the  Fcj's  of 

both  crystals  in  contact.  The  total  current 
undergoes  a  "refraction"  at  the  boundary.  In 
short,  boundaries  distort  the  current  flow 
and  the  distortion  extends  on  distances  of 
the  order  X  on  both  sides.  Another 
possibility  is  discussed  in  Ref.25:  the 
spontaneous  creation  of  vortex  chains  at 
sufficiently  long  boundaries  may  lead  to  an 
extra  dissipation  whenever  the  persistent 
current  crosses  the  interface. 

Some  support  for  the  presentation  of 
this  work  relates  to  Department  of  Navy 
Grant  N00014-89-J-1837  issued  by  the 
Office  of  Naval  Research.  The  US  Govern¬ 
ment  has  a  royalty-free  license  throughout 
the  world  in  all  copyrightable  materia) 
contained  herein. 
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In  the  absence  of  a  complete  microscopic  theory  of  the  high-temperature  copper-oxide  superconductors, 
several  phenomenological  theoretical  approaches  are  being  used  to  describe  the  intragranular  magnetic 
structure  in  the  mixed  state.  The  purpose  of  mis  paper  is  to  review  these  approaches  and  to  point  out  where 
further  developments  are  needed. 


1.  INTRODUCTION 

Since  the  discovery  of  the  high-temperature  oxide 
superconductors,1  these  materials  have  been  known 
to  have  magnetic  properties  characteristic  of  type-ll 
superconductors.  Although  the  properties  of  bulk 
sintered  samples  of  these  materials  are  greatly  com¬ 
plicated  by  the  effects  of  granularity,2  we  ignore  such 
effects  in  this  paper  and  focus  on  how  to  theoretically 
describe  the  intragranular  magnetic  structure  in  ideal 
untwinned  single  crystals,  which  are  known  to  be 
highly  anisotropic.  We  first  discuss  two  main  ap¬ 
proaches.  One  is  to  treat  the  magnetic  structure 
within  either  the  London  theory3  4  or  the  Ginzburg- 
Landau  theory5-'6  with  an  anisotropic  effective  mass 
tensor,  and  the  other  is  to  model  a  layered  supercon¬ 
ductor  as  a  stack  of  Josephson-coupled  supercon¬ 
ducting  blocks  or  layers.17'20  We  then  describe  the 
connections  between  these  two  approaches  and 
briefly  discuss  desired  extensions  of  the  theory. 

2.  ANISOTROPIC  EFFECTIVE-MASS  THEORY 

The  simplest  way  to  account  for  anisotropy  effects 

in  the  mixed  state  of  the  high-temperature  supercon¬ 
ductors  is  to  express  the  anisotropy  in  terms  of  a  phe¬ 


nomenological  effective-mass  tensor  in  the  framework 
of  the  London3'4  or  Ginzburg-Landau5'16  theories, 
which  normally  are  taken  to  be  isotropic.  In  the  refer¬ 
ence  frame  aligned  with  the  principal  axes,  this  mass 
tensor  is  diagonal,  and  the  diagonal  elements  m,  (/  = 

l ,  2,  3  =  a,  b,  c)  are  normalized13'14  such  that 

m,  m^m-i  =  1 .  The  penetration  depths  A,  =  Avm,  de¬ 
scribe  the  decay  of  components  of  the  supercurrent 
along  the  principal  directions,  and  the  corresponding 
coherence  lengths  4i  -  sA'm,  characterize  the  spatial 
variation  of  the  order  parameter  along  these  direc¬ 
tions.  The  scalars  A  =  (AiA2A3)1/3  and  4  =  (£is2s3)1/3 
enter  the  Ginzburg-Landau  expressions  for  the  bulk 
thermodynamic  critical  field  Hc  =  tp0Ai']2)2rc\$  (a 
scalar)  and  the  Ginzburg-Landau  parameter  xr=  A/£. 

Recent  vortex-lattice  decoration  experiments21'22 
have  shown  that  in  YBa2Cu307-s  the  penetration 
depths  are  approximately  in  the  ratios  Aa:At,Ac  = 

1 .2:1 :5.5.  These  results  thus  give  ma  -  0.4,  mf,~  0.3, 
and  mc  =  8.8.  To  explain  these  results  qualitatively, 
we  can  say  that  the  penetration  depths  Aa  and  A b  are 
much  smaller  than  Ac  because  the  screening  currents 
along  the  CuC>2  planes  (in  the  a  or  b  direction)  flow 
much  more  easily  than  in  the  c  direction.  Moreover, 
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the  Office  of  Naval  Research.  The  United  States  Government  has  a  royalty-free  license  throughout  the  world  in 
all  copyrightable  material  contained  herein. 


Ao  in  YBa2Cu307-5  is  somewhat  smaller  than  Aa  be¬ 
cause  current  parallel  to  the  CuC>2  planes  (lows  more 
easily  along  the  CuO  chains  (in  the  b  direction)  than 
perpendicular  to  these  chains  (in  the  a  direction). 

It  is  important  to  realize  that  the  penetration  depth 
depends  upon  the  direction  of  the  screening  super¬ 
current  density  j,  not  upon  the  direction  of  the  local 
magnetic  field  b.4-13-16  Suppose,  for  example,  an 
untwinned  single  crystal  is  cleaved  normal  to  a  prin¬ 
cipal  axis,  so  as  to  expose  a  flat  surface  in  the  yz- 
plane.  Application  of  a  weak  applied  magnetic  field 
Ha  along  a  principal  axis  parallel  to  the  surface  gen¬ 
erates  a  screening  current  density  j  along  the  other 
principal  axis  that  is  parallel  to  the  surface.  This  com¬ 
ponent  of  j  and  its  corresponding  local  magnetic  field 
b  {(4z/c)l  =  V  x  b]  decay  exponentially  as  exp(-x/A/), 
where  x  is  the  distance  into  the  sample  and  the  sub¬ 
script  i  refers  to  the  axis  along  which  j  flows. 

Within  Ginzburg-Landau  theory  with  an  anisotrop¬ 
ic  mass  tensor,5"16  the  magnetic  field  penetration  is 
governed  by  Maxwell's  equations  plus  the  second 
Ginzburg-Landau  equation, 

;,=  -(c/2/4  7rAy2)as/-.  (1) 

This  equation  gives  the  connection  between  the  /th  {/' 
=  1 , 2,  3  =  x,  y,  z)  component  of  the  supercurrent  den¬ 
sity  j  and  the  Ah  component  of  the  gauge-invariant 
vector  potential  as  =  a  +  ($o/2;t)Vy.  Here,  a  is  the 
vector  potential  (b  =  V  x  a)  and  yis  the  phase  of  the 
reduced  order  parameter  fexp(ry),  which  nan  be  re¬ 
garded  as  the  Cooper-pair  wave  function. 

For  a  vortex  centered  on  the  x-axis  (a-axis),  the 
London-model13’23  expression  for  the  magnetic  field, 
derived  assuming  f=*  1  outside  the  core,  is 

bx(y.  z)  -  (fc/2*Ai>Ac)Ko {[(yMc)2+(z/Ab)2]1/2).  (2) 

where  Ko  is  the  modified  Bessel  function  of  order 
zero.  The  decay  with  distance  along  the  y-axis  ( b- 
axis)  is  approximately  exponential,  with  decay  length 
Ac.  The  reason  that  the  decay  length  is  Ac  and  not  A  p 
is  that  the  decay  length  is  associated  with  the  direc¬ 
tion  of  the  screening  currents,  and  on  the  y-axis  these 


currents  point  in  the  z  direction.  Similarly,  along  the 
z- axis  (c-axis)  the  screening  currents  flow  parallel  to 
the  b  direction  and  hence  the  magnetic  field  decays 
exponentially  with  decay  length  A&.  Contours  of  con¬ 
stant  bx  and  streamlines  of  j  are  ellipses:  (y/Ac)2  + 
(zJ\b)2  =  const. 

The  core  of  the  London-model  vortex  is  defined  by 
the  ellipse  (y/Ac)Mz/Ab)2  =  (^Ac)2  =  (^A&)2  =  m,A2 
=Mk32.  Here  k=  A/£  is  the  Ginzburg-Landau  param¬ 
eter  and  *ra  =  K/^Jma  is  an  effective  nrfor  a  vortex 
aligned  along  the  a  direction.  The  semi-major  axis 
along  the  y-axis  (b-axis)  is  $b,  and  the  semi-minor 
axis  along  the  z-axis  (c-axis)  is  4c-  Inside  the  core.  bx 
=  const  and  f  =  0  according  to  the  London  model. 
Calculation  of  the  energy  per  unit  length  of  the  vortex 
along  the  a-axis,  £ia  =  (0o/4jr)Hc1a,  accounting  for  just 
the  magnetic  field  energy  density  tPi&rc  and  the  ki¬ 
netic  energy  density  {2nlc2){\b2jf+\c2jy1)  yields  the 
result,  valid  only  to  logarithmic  accuracy, 

Wcia  =  (tfto^rrApAcJInKa.  (3) 

A  more  accurate,  analytic  expression  for  the  mag¬ 
netic  field  of  an  isolated  vortex  in  the  vicinity  of  the 
core,  as  well  as  for  ci,  and  Hc\ ,,  can  be  obtained  by 
extending  the  variational  model  of  Ref.  24  to  the 
anisotropic  case.  In  this  approach,  one  takes  a  varia¬ 
tional  trial  function  for  the  order  parameter  of  the  form 
f=  p'/Ff',  where  p’=  [(y/Ac)2+(z/At>)2]1/2,  Ft'  =  (4v2 
4-(y/Ac)2+(z/At>)2]1/2,  and  /j/is  a  variational  core  ra¬ 
dius  parameter  of  order  1  Ik.  With  this  Ansatz,  the  so¬ 
lution  of  the  second  Ginzburg-Landau  equation 
is23"24 

bx(y,  z)  =  (<f*0/2^A£>Ac)Ko(H)  (4) 

when  k  »  1 .  When  this  result  and  f  *  p’/fl'  are  sub¬ 
stituted  back  into  the  Ginzburg-Landau  free  energy 
functional  to  minimize  the  energy  per  unit  length  of 
vortex,  all  the  integrals  can  be  performed  analytically 
and  the  energy  is  minimized  with  4v'  =  V 2/ira  if  *ra  » 
1,  yielding 


Ffcia=  (&/4/rA(>Ac)[Ko(V2/  Ka)  +  3/4), 


(5) 


or,  keeping  only  the  leading  terms  in  an  expansion  in 
powers  of  1  Ik3  , 

Hcia=  (0o/4^c)(lmr3+ 0.52).  (6) 

The  above  variational  calculation,  however,  slightly 
overestimates  the  energy  per  unit  length  of  the  vortex. 
As  shown  in  Ref.  25,  numerical  solutions  of  the 
Ginzburg-Landau  equations  show  that  the  constant 
term  is  not  0.52  but  0.497.  Note  that,  if  the  vortex  is 
centered  along  the  y-  or  z-axis,  the  corresponding  ex¬ 
pressions  for  b,and  Hc i,can  be  obtained  by  cyclic 
permutation  (x— >y->z-*x  or  a->b— >c-+a). 

For  a  discussion  of  the  more  complicated  case  for 
which  the  vortex  is  not  along  any  of  the  principal  axes, 
the  reader  is  referred  to  Ref.  26. 

At  the  upper  critical  field,  the  linearized  Ginzburg- 
Landau  equations  can  be  solved  for  arbitrary  angles 
(0,<p)  of  the  applied  field  H  =  H  ($  sinOcostp  *  P  sin# 
sind  +  t  cose).  The  magnitude  of  the  upper  critical 
field  (at  which  bulk  superconductivity  is  quenched)  is9 
Hc2(0,<p)  =  21'2kwK0.*)Hc«  (do/2^7)(me,K#.d)]-1/2. 
where  <eii{S.<P)  =  *trne//< 0,d)]- 1/2  and 

meii{Q,<p)  =  ma  sin20cos20+  mt,  sin20sin2d 

+  mccos20.  (7) 

When  the  applied  field  is  directed  along  the  a-,  b-, 
and  c-axes,  we  have  Hcza  =  frJ2itSbiSc.  HcZb  B 
do/2 ttScSa,  and  Hc2c  =  respectively,  as  can 

be  shown  with  the  help  of  mi  012073  =  =  f  and 

Si  =  &M- 

3.  JOSEPHSON-COUPLED  BLOCKS  OR  LAVERS 
Various  consequences  of  modeling  the  high-tem¬ 
perature  superconductors  as  Josephson-coupled 
blocks,  layers,  or  subgrains1 7-20  are  currently  under 
theoretical  and  experimental  investigation.  In  order  to 
bring  out  how  magnetic  fields  penetrate  a  supercon¬ 
ductor  in  such  a  picture,  we  analyze  here  a  general 
model,  sketched  in  Fig.  1,  consisting  of  a  periodic  ar¬ 
ray  of  Josephson-coupled  blocks.  The  blocks,  rect¬ 
angular  parallelepipeds  of  dimensions  ai\  a2'.  and 
aj'  along  the  x,  y,  and  z  (a,  b,  and  c)  directions,  are 


arranged  on  an  orthorhomoic  lattice  with  lattice  pa¬ 
rameters  ai.  32,  and  a 3.  The  areas  of  the  junctions 
between  the  blocks  normal  to  i.  P,  and  i  are  Ai'  = 
a2'a3',  A-y  =  a3'ai',  and  A3'  =  ai'a2\  while  the  corre¬ 
sponding  areas  of  the  sides  of  the  unit  cell  are  At  = 
a2a3,  A2  =  a3ai ,  and  A3  =  ai  a2.  The  blocks  are 
aligned  with  their  symmetry  axes  along  the  x,  y,  and  z 
directions,  such  that  the  corresponding  intrinsic  pen¬ 
etration  depths  are  A10 (T),  A2o(7),  and  A3o(7),  while 
the  intrinsic  coherence  distances  are  £io(7),  Szo(T). 
and  SzoCH-  We  assume  that  a,  '  «  A«,  such  that  the 
local  magnetic  field  within  the  unit  cell  penetrates 
each  block  almost  uniformly.  The  blocks  are  Joseph¬ 
son-coupled  with  coupling  energy  Ej/(T)  =  (tr/2e)ljo{T) 
across  the  junctions  normal  to  the  /  direction.  The  su¬ 
percurrent  in  the  /  direction  across  such  a  junction  is  /, 
=  /,osind)t,  where  !,o (7)  is  the  maximum  dc  josephson 
current,27 

lo(T)  =  (xA(T)/2eR,n]tanh(A(7)/2ke71,  (8) 

A(T)  is  the  temperature-dependent  gap  parameter,  R,n 
is  the  normal-state  tunneling  resistance  across  the 
junction,  and  Ay;  is  the  gauge-invariant  phase  differ¬ 
ence  across  this  junction. 

Regardless  of  the  strength  of  the  Josephson  cou¬ 
pling,  the  local  magnetic  field  varies  slowly  on  the 
scale  of  a,.  We  wish  to  obtain  a  differential  equation 
describing  how  superconductivity  affects  the  slow 
spatial  variation  of  B,  the  magnetic  flux  density  aver¬ 
aged  over  the  unit  cell.  This  equation,  in  combination 
with  Ampere's  law  (4n/c)J  =  V  x  B.  where  J  is  the 
current  density  averaged  over  the  unit  cell,  will  yield 
useful  expressions  for  the  supercurrent  penetration 
depths  A,  along  the  three  principal  directions.  Using  a 
procedure  similar  to  that  used  in  Ref.  28.  we  obtain 
the  desired  equation  by  considering  the  magnetic  flux 

S3 (x,y)A3  in  the  z  direction  through  the  rectangle  of 
area  A3  =  aia2  centered  at  P(x,y).  This  flux  can  be 
expressed  as  a  line  integral  of  the  vector  potential  a 
around  the  rectangular  contour  1-2-3-4-5-6-7-8-9-10- 
1 1  -1 2-1  shown  in  Fig.  1 .  Within  the  blocks  we  first  re¬ 
place  a  by  as-  (&>/ 2x)Vy.  We  next  make  use  of 


I 


FIGURE  1 

Josephson-coupled-block  model. 


•  (4;rA io2lc)ji.  where  yV  =  J,Ah  and  the  magnetic  field 
is  assumed  to  be  sufficiently  weak  that  f-  1.  Ex¬ 
pressing  differences  of  J,  in  adjacent  blocks  in  terms 
of  partial  derivatives,  we  then  obtain  the  following 
fluxoid  quantization  condition, 

S3  =  '  (0o/2;r-43)(/}}s  4  +  4/3,7  +  Aft,  10  +  4712,1) 

+  (4nlc)[(kw2a\A\la\A\)dJt(x,y)ldy 

*  (^ir^azA^azA^dJzix.yydx],  (9) 

where,  for  example,  the  gauge-invariant  phase  differ¬ 
ence  across  junction  9, 1 0  is 

10 

Aft.io  =  n  -  no  ■  (2rfM  fdi-a.  (10) 

9 

the  line  integral  being  taken  directly  across  the  junc¬ 
tion. 

Since  473,4  *  -Ay\{x,y+azl2),  476,7  = 

-472(x-ai/2,y),  473,10  -  47i(x,y-a2/2),  and  4712,1  = 
472(x+ai/2,y),  we  also  may  write 

S3  -  (£y'2/&43)[a2c)47'i(x,y)/dy-  ai9472(x,y)] 

+(4  n/c)[(  A 1  o2ai  'A \!a\  A 1 ')d  J\  (x,y)/3y 

-  (A.202a2'42/a242')cU2(x,y)/8x].  (11) 


For  weak  fields  and  small  current  densities,  we  may 
linearize  the  Josephson  current  equation  (replace 
sin4yby  4 7).  Substitution  of  J,(x,y)  =  J/oAyiix.y), 
where  Jo(r,>  =  /<}( 7V4,.  into  the  above  equation  then 
yields 

S3  =  (4a/c)(Ai20Ji/0y  -  X22dJzldx).  (1 2) 

where 

A/2  =  Xo2aj‘A/ajAi’  +  cipo/8i^a^Jio.  (f3) 

Corresponding  expressions  for  Si  "md  can  be 
obtained  from  the  above  by  cyclic  permutation.  From 
these  expressions  and  Ampere’s  law,  one  may  show 
that  if  J  flows  along  only  one  of  the  three  principal 
axes,  only  the  A,  along  that  axis  determines  its  spatial 
variation.  For  example,  if  J  =  Jjji ,  then 

J2jd  xX-  +  =  J3/A32.  (14) 

To  discuss  the  limiting  cases  of  the  expression  tor 
A,,  note  that  if  the  junction  barrier  thicknesses  shrink  to 
zero  and  then  and  A/-+A,.  Thus 

A,-» A,o  in  the  strong-Josephson-coupted  limit.  On  the 
other  hand,  if  J/o  «  cipoA/IQrfiX^Aai',  we  have  A/  = 
(c0O/'842a,Jio)1/2  in  the  weak-Josephson-coupling 
limit. 

If  the  intrinsic  normal-state  resistivity  within  a  given 
block  along  the  /-axis  is  pa.  then  the  effective  resistiv¬ 
ity  of  the  array  of  blocks,  which  are  connected  in  se¬ 
ries  with  the  tunnel  resistances  Rm.  is 

pi  =  poaj'A/aiAj’  ♦  RinA/aj.  ( 1 5) 

In  the  limit  when  Rm  «  pnai'/Aj',  which  can  occur  in 
the  strong-Josephson-coupled  limit,  the  junction  re¬ 
sistance  makes  a  small  contribution  to  the  overall  re¬ 
sistivity,  while  in  the  opposite  limit  when  Rm  » 
poai'/Aj',  the  junction  resistance  dominates  the  ex¬ 
pression  for  pi. 

If  the  superconducting  blocks  are  in  the  clean  limit, 
there  is  no  relation  between  ratios  of  the  Ao’s  and  the 
pa's.  Similarly,  for  an  array  of  clean-limit  blocks, 


there  is  no  relation  between  ratios  of  the  A/s  and  the 
p/s.  However,  in  the  weak-coupling  limits  for  which 
Rm  »  poa/'/A'and  Jq  «  ap0Al'J8,'fiX/o2Alal',  the  Am- 
tegaokar-Baratoff  theory27  gives  the  following  rela¬ 
tion  between  X,  and  p, : 

A,2  =  frc2p/4rt24(T)tanh[4(T)/2/cen.  (16) 

which  agrees  exactly  with  the  corresponding  result 
obtained  from  the  BCS  theory  in  the  dirty  limit.29  30 
These  results  thus  explain  how  it  is  possible  for  a 
highly  anisotropic  high-temperature  superconductor 
to  behave  as  if  it  were  in  the  clean  limit  along  one  di¬ 
rection  and  in  the  dirty  limit  along  another. 

If  we  consider  the  case  for  which  the  blocks  are 
strongly  Josephson-coupled  along  the  a  and  b  direc¬ 
tions  but  weakly  coupled  along  the  c  direction,  we 
obtain  results  similar  to  that  of  the  closely  related 
lawrence-Doniach  model.’7  In  this  case  the  corre¬ 
sponding  penetration  depths  along  the  a,  b,  and  c  di¬ 
rections  are  given  by 


^  a2  =  A  ao2  33/33', 

(17a) 

Xb2  =  Xbo2azla3, 

(17b) 

Ac2  =  Aco233'/a3  +  ctpo/Bn^ajJzo, 

(17c) 

where  is  the  lattice  parameter  along  the  c  direction, 
a3'  is  the  effective  thickness  of  the  superconducting 
layer,  and  J30  is  the  maximum  Josephson  supercur¬ 
rent  density  perpendicular  to  the  layers. 

The  above  derivations  show  how  an  anisotropic 
penetration  depth  tensor  can  be  derived  for  an  ori¬ 
ented  array  of  Josephson-coupled  blocks  or  layers 
whose  small  dimensions  are  less  than  the  intrinsic 
penetration  depth.  The  magnetic  structure  in  such  an 
array  strongly  resembles  that  in  a  homogeneous  su¬ 
perconductor  described  by  the  London  or  Ginzburg- 
Landau  theories  with  an  anisotropic  effective  mass 
tensor.  Far  from  the  vortex  cores  the  behavior  is 
nearly  the  same  in  the  two  descriptions,  but  differ¬ 
ences  may  occur  in  the  vicinity  of  the  vortex  "core."  In 
a  homogeneous  superconductor  the  1/r  dependence 
of  the  superfluid  velocity  around  the  vortex  axis  is  cut 
off  by  order-parameter  suppression  (/<*  r)  in  the  core, 


such  that  the  current  density  becomes  proportional  to 
r  close  to  the  vortex  axis.  Whether  or  not  strong  order- 
parameter  suppression  in  the  core  occurs  in  a 
Josephson-coupled  array  depends  upon  the  size  of 
the  effective  coherence  distance  relative  to  a/.  The 
results  of  Ref.  30  suggest  that,  while  may  be  close 
to  si/  o  with  sufficiently  strong  Josephson  coupling,  the 
appropriate  expression  in  the  weak-Josephson-cou- 
pting  limit  is  given  near  Tcby 

= ’t1rAHT)UeZkBTcH<2(T)pi.  (18) 

The  size  of  the  vortex  "core"  is  then  given  by  the 
larger  of  ^  or  all.  If  <  a/2,  the  suppression  of  the 
order  parameter  of  the  block  is  very  small,  the  vortex 
axis  sits  in  the  Josephson  junction  between 
blocks,18'19  and  the  lower  cutoff  radius  (in  expres¬ 
sions  such  as  that  for  Hc i)  is  all,  not  & . 

4.  CONCLUSION 

In  this  paper,  I  have  discussed  several  ap¬ 
proaches  being  used  to  describe  the  intragranular 
magnetic  structure  in  the  mixed  state  of  the  high-tem¬ 
perature  superconductors.  Further  extensions  of 
these  models  are  needed  to  provide  a  clear  picture  of 
the  crossover  between  the  regime  in  which  the  mag¬ 
netic  structure  varies  on  length  scales  A,  and  ^  that 
are  both  larger  than  the  linear  dimensions  of  the  dis¬ 
crete  building  blocks  of  the  superconductor,  and  the 
regime  in  which  the  smallest  length  scales  of  the 
magnetic  structure  are  set  by  the  sizes  of  these 
building  blocks. 
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Recent  photoemission  studies  of  Y  123  and  81  2212  performed  with  high  energy  and  angular 
resolution  have  provided  detailed  Information  on  the  nature  of  the  states  near  the  Fermi 
level.  Measurements  of  the  superconducting  gap,  band  dispersion,  and  the  density  of 
states  near  the  Fermi  level  In  the  normal  state  all  support  a  Fermi  liquid  description  of 
these  materials. 


1.  INTRODUCTION 

Historically,  photoemission  has  had  alxed 
success  in  determining  the  electronic  structure 
of  the  high  temperature  superconductors.  More 
recently,  it  has  been  shown1  that  It  Is 
important  to  cleave  YBajCu,0T_4  (Y  123)  samples 
at  low  temperature  to  ensure  that  the  surface 
probed  by  photoemission  Is  characteristic  of 
the  bulk.  The  BijSrjCaCu.,0,  (B1  2212)  class  of 
superconductors  Is  much  more  stable  at  elevated 
temperature,  but  reduced  teaperatures  prolong 
the  useful  life  of  a  surface.  Additionally, 
reduced  Fermi  level  broadening  is  important  In 
studying  the  details  of  band  dispersion  near 
the  Fermi  level.  Spectra  can  be  further  simp¬ 
lified  by  angular  (k  space)  resolution,  high 
kinetic  energy  resolution,  or  the  application 
of  all  three  at  the  same  time. 


2.  EXPERIMENTAL 

Measurements  were  done  at  the  Synchrotron 
Radiation  Center  ,  Stoughton,  Wisconsin  on  the 
Aaes/Hontana  State-ERG/Seya  beamline2.  The 
photoemission  spectrometer  has  two  degree 
angular  resolution  and  Is  capable  of  20  meV 
electron  resolution.  The  combined  system 
resolution  In  Figs.  3  and  4  was  32  aeV. 

3.  RESULTS 

The  spectra  In  Fig.  1  are  taken  In  normal 
emission.  In  this  geometry  a  shift  In  the 
binding  energy  of  a  peak  as  a  function  of 
photon  energy  represents  dispersion  along  the 
direction  normal  to  the  ab  plane.  As  can  be 
seen  In  Fig.  la  for  B1  2212,  there  Is  minimal 
dispersion.  The  primary  effects  are  a  small 
shift  of  the  3.5  eV  peak  to  lower  binding  energy 
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FIGURE  1 

Photoemlsslon  energy  distribution  curves  taken 
in  normal  emission  for  Bi  2212  (a)  and  Y  123  (b). 
The  samples  were  cleaved  and  maintained  at  20K. 
The  spectra  are  normalized  to  incident  photon 
flux. 

with  increasing  photon  energy  and  a  larger  shift 
of  the  4.8  ev  peak  to  higher  binding  energy. 

The  small  amount  of  dispersion  in  the  z  direc¬ 
tion  is  not  surprising  in  this  highly  two  dimen¬ 
sional  material.  Dispersion  in  less  well 
defined  peaks  is  masked  by  the  very  strong 
intensity  modulations  due  to  matrix  element 
effects.  The  enhancement  near  the  Fermi  level 
and  at  1.5  ev  for  photon  energies  near  18  eV  are 
only  one  of  the  strong  selection  rule  effects. 

Fig.  lb  is  a  similar  set  of  spectra  for 
Y  123.  Again,  the  two  dimensionality  of  this 
material  results  in  minimal  dispersion  in  this 
geometry.  There  are  indications  of  dispersion 
lower  in  the  valence  band,  but  the  dominant 
features  are  the  very  strong  matrix  element 
effects.  From  these  spectra,  and  others  taken 
at  other  emission  angles,  it  is  clear  that  the 
strong  enhancement  of  the  1.2S  ev  peak  at  a 
22  ev  photon  energy  is  a  matrix  element  effect 
and  not  an  oxygen  resonance. 


Dispersion  is  much  larger  along  directions 
in  the  ab  plane.  Fig.  2  are  spectra  for 
different  emission  angles  and  hence  different  k 
values  along  approxlmatly  the  r-S  direction.  An 
initial  state  at  1.25  eV  splits  with  one  branch 
remaining  at  about  1.25  ev  and  the  other  moving 
deeper  away  from  the  zone  center.  Additional 
dispersion  is  seen  in  the  lower  part  of  the 
valence  band.  Finally,  a  band  is  seen  disper¬ 
sing  through  the  Fermi  level.  In  both  Fig.  1 
and  2,  because  of  the  moderate  resolution 
(150  asV),  esissior.  near  the  rermi  level  is 
largest  when  the  band  is  actually  well  below  E^. 
High  resolution  data  taken  at  the  same  time  as 
the  spectra  In  Fig.  2  show  that  the  Fermi  level 
crossing  is  at  12  degrees. 
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FIGURE  2 

Angle  resolved  photoemission  spectra  along 
approximately  the  T-S  direction,  normalized  to 
incident  photon  flux.  The  angle  is  relative  to 
the  sample  normal. 

For  Bi  2212,  emission  near  the  Fermi  level 
can  be  seen  in  much  greater  detail  in  the  high 
resolution  spectra  of  Fig.  3.  As  predicted  by 
mil  band  structure  calculations,  there  is 
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FIGURE  3 

High  resolution  angle  resolved  spectra  along 
approximately  the  T-H  direction,  normalized  to 
incident  flux. 

minimal  emission  from  the  Fermi  level  near  the 
zone  center.  The  structure  that  Is  seen  near 
normal  emission  in  the  figure  probably  comes 
from  randomly  oriented  damaged  areas  at  the 
edges  of  the  sample.  Increasing  angles  sample 
a  line  near  the  zone  diagonal  (T-H).  A  band  Is 
seen  dispersing  towards  the  Fermi  level  from  at 
least  300  aeV  below  E^  and  then  moving  away 
from  the  Fermi  level.  This  follows  the 
calculated  band3  qualitatlvly  but  the  measured 
band  appears  ouch  closer  to  the  Feral  level 
than  predicted.  This  Is  In  qualitative 
agreement  with  Takahashi  et.al. 

The  measurements  shown  In  Fig.  3  were  taken 
at  20K,  and  the  shape  of  the  band  near  the  Fermi 
level  is  distorted  by  the  existence  of  a 
superconducting  gap.  Data  reported  In  Ref.  S 
were  taken  along  a  line  parallel  to  that  of 
Fig.  3,  and  Include  a  point  which  Is  one  of  the 
few  places  where  a  band  in  the  normal  state 
crosses  the  Fermi  level.  The  spectra*  In  the 
normal  and  superconducting  state  are  shown  In 


the  upper  part  of  Fig.  4.  rtodellng  these 
spectra  gave  2&/kTc  *7.  twice  the  weak  coupling 
BCS  value.  In  the  bottom  half  of  Fig.  4,  the 
same  normal  state  spectrum  Is  compared  to  a 
platinum  foil  on  the  sample  holder.  The 
instrumental  resolution  (32  meV)  is  comparable 
to  the  Fermi  level  width  at  this  temperature 
(34  meV).  The  measured  Fermi  levels  are 
equivalent  to  within  the  resolution  and  noise 
of  the  experiment.  A  report  by  these  authors 
of  a  more  extensive  experiment  on  the  normal 
state  Feral  level  properties  Is  in  preparation, 
with  the  goal  of  placing  quantitative  limits  on 
the  possible  value  of  RVB  parameters. 
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FIGURE  4 

Upper:  High  resolution  spectra  showing 
modifications  to  the  density  at  the  Fermi  level 
as  a  result  of  a  superconducting  gap.  The 
solid  curve  Is  for  20K,  and  the  dotted  curve  Is 
for  90K;  above  the  T  of  82K. 

Lower:  The  same  norial  state  density  compared  to 
a  Pt  density  of  states  (solid  line)  at  90K. 

4.  CONCLUSIONS 

High  energy  and  angular  resolution  combined 
with  low  temperature  is  giving  an  increasingly 
detailed  picture  of  the  high  temperature 
superconductors.  Some  of  the  detail,  such  as 
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The  magnetic  properties  of  the  half-filled  one  band  Hubbard  model  is  studied  within  the  random 
phase  approximation  (RPA)  above  the  mean  field  spin-density-wave  (SDW)  background.  The  spin 
wave  spectrum  and  the  sublattice  magnetization  are  shown  to  interpolate  smoothly  between  the  weak 
and  strong  coupling  limits.  A  hole  in  the  SDW  background  is  self-trapped  in  a  region  of  reduced 
SDW  amplitude  and  forms  an  object  which  we  call  the  spin  bag.  The  attraction  among  the  two  hoies 
results  due  to  the  tendency  to  share  the  region  of  reduced  SDW  amplitude.  A  detailed  calculation 
based  on  the  RPA  above  the  SDW  background  support  this  picture.  The  symmetries  and  the  nodes 
of  the  resulting  pairing  gap  are  discussed. 


1.  INTRODUCTION 

One  of  the  characteristic  features  of  the  high  Tc 
superconductors  is  the  antiferromagnetic  correlations 
present  both  in  the  insulating  and  in  the  metallic 
phases.1  Many  theories  are  thus  based  on  the  point 
of  view  that  this  observed  antiferromagnetic  correla¬ 
tion  is  intimately  related  to  the  occurance  of  the  su¬ 
perconductivity.  There  are  two  different  ways  to  ap¬ 
proach  the  magnetic  phase  based  on  the  single  band 
Hubbard  model,  the  localized  Mott  limit  where  the 
problem  can  be  mapped  onto  a  Heisenberg  model  | 
and  the  itinerant  limit  where  the  magnetic  order  ! 
arises  due  to  the  fermi  surface  nesting.  We  shail 
present  a  quantitative  analysis  based  on  the  RPA 
above  the  SDW  background  to  show  that  physi¬ 
cal  properties  actually  interpolate  smoothly  between 
these  two  seemingly  opposite  limits.2  The  doped 
metallic  phase  can  also  be  approached  from  two  dif¬ 
ferent  directions,  with  one  starting  from  the  para¬ 


magnetic  phase  where  one  postulates  a  new  spin  liq¬ 
uid  like  magnetic  order  and  the  other  starting  from 
the  antiferromagnetic  phase  where  one  assumes  that 
the  length  scale  of  the  antiferromagnetic  order  to 
be  larger  than  the  other  length  scales  in  the  prob¬ 
lem.  One  expects  that  physical  properties  should 
also  interpolate  smoothly  between  these  two  limits 
although  a  quantitative  demonstration  is  still  lack¬ 
ing.  In  our  approach,2  we  shall  discuss  the  problem 
starting  from  the  antiferromagnetic  phase. 

2.  MAGNETIC  PROPERTIES  AT  HALF-FILLING 

At  half-filling  and  strong  coupling  ( U  >  t),  the 
Hubbard  model  maps  onto  the  Heisenberg  model 
with  J  ~  t2/U.  The  numerical  investigations  show 
that  the  ground  state  is  N6el  ordered.3  The  spin 
wave  velocity  is  of  the  order  of  J  and  the  quan¬ 
tum  fluctuations  reduce  the  sublattice  magnetiza- 
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tion  to  0.6)j.g  from  the  classical  value  of  l^tg.4 
the  weak  coupling  limit  ( U  <§:  t ),  the  fermi  surface 
is  a  perfect  square  and  this  nesting  property  leads 
to  an  instability  of  the  free  fermi  surface  and  to  the 
formation  of  a  SDVV  ground  state.  In  this  limit, 
the  spin  wave  velocity  is  of  the  order  of  the  fermi 

i 

velocity  t  and  the  sublattice  magnetization  is  expo-  j 
nentially  small.  Clearly  the  ground  states  in  both 
limits  have  the  same  symmetry  and  other  qualita¬ 
tive  features,  the  question  is  how  physical  quanti¬ 
ties  interpolate  between  these  two  limits.  We  per¬ 
formed  a  calculation  based  on  the  RPA  above  the 
mean  field  SDW’  ground  state  to  extract  the  spin 
wave  velocity  from  the  gapless  pole  in  the  trans¬ 
verse  susceptibility."  The  resulting  spin  wave  veloc¬ 
ity  indeed  interpolates  smoothly  between  the  weak 
and  the  strong  coupling  limits.  Furthermore,  both 
the  mean  field  and  the  fluctuation  reduced  sublattice 
magnetization  can  be  worked  out  within  this  formal¬ 
ism  and  again  we  found  smooth  behavior  between 
the  two  limits.  Actually  the  fluctuation  reduced  sub¬ 
lattice  magnetization  we  obtained  compare  well  with 
the  Monte  Carlo  results  of  Hirsrh  and  Tang.3 

3.  THE  SPIN  BAG  AND  SUPERCONDUCTIV¬ 
ITY 

When  a  hole  is  injected  to  the  rigid  SDW  back¬ 
ground,  it  costs  a  minimal  energy  of  the  SDW  gap 
Asdw-  The  system  therefore  has  the  tendency  of  re¬ 
ducing  the  SDW  amplitude  in  the  vicinity  of  the  hole 
as  to  lower  the  self-energy.  The  gain  in  self-energy 
of  the  hole  is  balanced  by  the  cost  of  the  exchange 
energy  and  these  two  competing  effects  lead  to  the 
self-consistent  formation  of  an  object  which  we  call 
a  spin  bag.2  The  size  of  the  spin  bag  is  of  the  order 
of  the  SDW  coherence  length  £sDW  */&SDW 
is  about  three  to  four  lattice  spacing.  The  shape  of 
the  spin  bag  can  be  round  or  cigar  like  depending 
on  the  location  in  fermi  surface.  All  these  quali¬ 


tative  features  are  observed  In  the  numerical  varia¬ 
tional  calculation  by  Su.5  From  this  reasoning  one 
concludes  that  two  holes  attract  each  other  because 
of  the  tendency  to  share  the  region  of  reduced  SDW 
amplitude  as  to  lower  the  cost  of  exchange  energy. 
We  believe  that  this  is  the  basic  mechanism  by  which 
the  superconducting  pairing  occurs. 

The  above  qualitative  picture  is  supported  by  a 
detailed  dynamic  calculation  of  the  pairing  interac¬ 
tion  mediated  by  the  collective  modes  of  the  SDW 
background.2,6  These  collective  modes  are:  the  charge 
fluctuation,  the  spin  orientation  and  the  spin  ampli¬ 
tude  fluctuations.  There  is  no  significant  enhance¬ 
ment  in  the  charge  fluctuation  channel  since  one  is 
dealing  with  the  positive  U  Hubbard  model.  The 
spin  orientation  fluctuation  (the  spin  wave)  is  an  im¬ 
portant  low  energy  mode.  However  it  couples  to  the 
holes  with  a  vertex  which  vanishes  in  the  small  dop¬ 
ing  limit  and  can  be  neglected  in  the  Born  approx¬ 
imation.  The  most  important  mode  is  the  spin  am¬ 
plitude  fluctuation  which  couples  to  the  holes  with 
unit  strength  even  in  the  small  doping  limit  and  ac¬ 
cording  to  our  previous  argument,  it  is  responsible 
for  the  formation  of  the  spin  bag. 

The  qualitative  features  of  the  spin  amplitude 
fluctuation  are:  I)  attractive  for  small  momentum 
transfer  and  repulsive  for  momentum  transfer  near 
Q  =  (?r,ir)  the  nesting  wave  vector;  2)  the  spectral 
weight  is  peaked  around  Ag^vv.  which  serves  as  a 
frequency  cut-off  similar  to  the  Debye  frequency  in 
the  phonon  problem.  For  small  doping,  holes  form 
“pocket"  at  the  magnetic  zone  boundary,  concen¬ 
trated  either  around  (0, ir)  or  ( j,  §).  These  features 
of  the  hole  fermi  surface  combined  with  the  proper¬ 
ties  of  the  spin  amplitude  fluctuation  mode  lead  to 
a  pairing  gap  that  has  either  a  d- wave  or  p- wave  like 
symmetry  but  the  pairing  gap  has  no  nodes  at  the 
hole  fermi  surface  in  either  case. 


4.  CONCLUSION 


Starting  from  the  mean  field  SDVV  ground  state, 
we  demonstrated  that  RPA  fluctuations  describe  the 
magnetic  properties  correctly  in  both  the  localized 
and  the  itinerant  limits.  Upon  doping,  the  spin  am¬ 
plitude  fluctuation  is  the  dominant  mode  that  cou¬ 
ples  to  the  hole  and  leads  to  the  formation  of  the 
spin  bag.  The  properties  of  the  hole  fermi  surface 
and  the  spin  amplitude  mode  are  such  that  the  re¬ 
sulting  pairing  gap  is  nodeless  over  the  hole  fermi 
surface. 


ACKNOWLEDGEMENT 

We  would  like  to  acknowledge  helpful  discussions 
with  H.  Nakanishi,  D.J.  Scalapino  and  K.  Schonham- 
mer.  This  work  was  supported  in  part  (J.R.S.  and 
S.C.Z.)  by  the  NSF  grant  DMR85-17276  and  also 
(X.G.W.)  PHY82-17853,  supplemented  by  funds  from 
NASA. 

REFERENCES 

1.  G.  Shirane,  et  al.,  Phys.  Rev.  Lett.  59  (1987) 
1613;  J.  Transquada,  et  al .,  Phys.  Rev.  Lett.  60 
(1988)  156;  K.B.  Lyons,  et  al.,  Phys.  Rev.  B37 
(1988)  2353;  K.B.  Lyons,  et  al.,  Phys.  Rev.  Lett. 
60  (1988)  732;  R.J.  Birgeneau,  et  al.,  Phy3.  Rev. 
B38  (1988)  6614. 

2.  J.R.  Schrieffer,  X.G.  Wen  and  S.C.  Zhang,  Phys. 
Rev.  Lett.  60  (1988)  944;  Phys.  Rev.  B39 
(1989)  11663. 

3.  J.D.  Reger  and  A.P.  Young,  Phys.  Rev.  B37 
(1988)  5978;  J.E.  Hirsch  and  S.  Tang,  preprint. 

4.  P.W.  Anderson,  Phys.  Rev.  86  (1952)  694. 

5.  W.P.  Su,  Phys.  Rev.  B37  (1988)  9904. 

6.  Z.Y.  Weng,  T.K.  Lee  and  C.S.  Ting,  Phys.  Rev. 
B38  (1988)  6561;  G.  Vignale  and  K.S.  Singwi, 
Phys.  Rev.  B39  (1989)  2956. 


Probing  Magnetism  and  Superconductivity  of  High*rc  Systems  with  Positive  Muons 


Y.J.  Uemura’ ,  L.P.  Lel,  G.M.  Luke1,  B.J.  Sternlieb1, 

J.H.  Brewer2,  R.  Kadono2  R.F.  Kiefl2,  S.R.  Kreitzman2,  and  T.M.  Riseman2 

1.  Department  cf  Physics,  Columbia  University,  New  York  City,  New  York  100S 7,  U.S-A. 

1.  TRIUMF  and  Department  of  Physics,  University  of  British  Columbia,  Vancouver,  B.C.,  V6T 
2 A3,  Canada 

Abstract 

Recent  results  from  our  /iSR  experiments  on  high-Tc  systems  are  presented.  We  describe  studies 
of  static  magnetic  order  in  hole-doped  and  electron-doped  systems,  correlations  between  Tc  and 
n,jm"  (carrier  density  /  effective  mass)  in  hole-doped  superconducting  cuprates,  the  temperature 
dependence  of  the  penetration  depth  in  c-axis  oriented  sintered  ceramic  of  Y  Ba-iCu^Ot,  as  well  as 
the  detection  of  flux  pinning  phenomena  by  pSR. 


Introduction 

Positive  muon  spin  relaxation  (pSR)  has  been  very  use¬ 
ful  for  the  study  of  high- Tc  systems1-3.  With  this  technique, 
one  can  measure  static  magnetic  order,  even  with  very  small 
and  random  moments,  and  study  details  of  magnetic  phase 
diagrams.  The  magnetic  field  penetration  depth  in  supercon¬ 
ductors,  determined  directly  by  pSR,  has  reveaded  various 
aspects  of  superconducting  energy  gap,  carrier  density,  as 
well  as  flux  pinning.  Ln  this  paper,  we  review  recent  progress 
of  our  pSR  experiments  at  TRIUMF  (Vancouver)  to  demon¬ 
strate  these  features  of  pSR. 

Experimental  aspects 

The  principle  of  /jSR  technique  has  been  described  in 
earlier  review  papers1,2  and  other  literatures'*.  The  polar¬ 
ized  positive  muon  beam  is  stopped  in  a  specimen  (typically 
1  ~  2  cm  in  diameter,  1  ~  2  mm  thick),  and  time  his¬ 
tograms  of  muon  decay  events  are  recorded,  one  by  one,  by 
positron  counters  placed  at  forward  (F)  and  backward  (B)  di¬ 
rections  with  respect  to  the  initial  muon  spin  direction.  Since 
positrons  are  emitted  preferentially  along  the  muon  spin  di¬ 
rection,  the  time  evolution  of  the  muon  spin  direction  and 
polarization  can  be  directly  observed  via  the  F/B  angular 
asymmetry  of  the  positron  time  spectra.  This  measurement 
can  be  performed  either  in  zero  field  (ZF-pSR)  or  by  applying 
a  transverse  external  magnetic  field  (TF-pSR).  After  stop¬ 
ping  almost  instantaneously  (typically  within  ~  10-11  sec) 
at  an  interstitial  site,  the  generally  stays  at  rest  in  oxides 
during  its  life  (mean  life  time  of  is  2.2  psec),  probing  mag¬ 
nitude.  randomness,  and  fluctuations  of  local  magnetic  fields 
at  the  interstitial  site,  analoguously  to  an  implanted  proton. 

Magnetic  order  —  hole-doped  systems 

We  previously  reported  4SR  studies  in  LaiCuOt-^  (ref. 
5)  and  VflujCujOx  (ref.  6)  which  established:  (1)  The 
Neel  temperature  is  reduced  sharply  with  increasing  Oxy¬ 
gen  concentration,  and  static  mgnetic  order  disappears  near 
the  concentration  where  the  superconductivity  appears;  (2) 
the  static  moment  of  Cu  in  the  antiferromagnetic  phase  does 


not  depend  much  on  the  Neel  temperature.  There  were,  how¬ 
ever,  a  report  by  VVeidinger  et  al.7  that  the  static  magnetic 
order  remains  deep  in  the  suprconducting  phase,  e  g.,  even  in 
La  1. gs T r0  ! 5 CuCV  To  clarify  whether  or  not  this  is  the  in¬ 
trinsic  propertiy,  we  have  performed8  zero-field  /rSR  measure¬ 
ments  on  a  Lai  gsSro  isCuOi  specimen  made  by  spray-dry 
method  as  well  as  on  carefully  annealed  sintered  specimens 
of  Y BaiCu3Ox  at  T  =  0  02  K. 

Zero  Field.  T  =  20  mK 


Fig.  /.  Zero-field  pSR  spectra  at  T  =  20  mK  in  Y  Ba^C U3O 
with  x  =  0.27  (non- superconducting),  x  =  6.54  (super con¬ 
ducting),  and  in  Lai.asSro .isCuOi  (from  Kiefl  et  elf). 


Figure  1  shows  the  time  dependence  of  the  muon  spin 
polarization  in  these  systems.  In  non-superconducting 
Y Ba^Cu^Oe.n,  the  muon  spin  precesses  around  the  static  in¬ 
ternal  field  from  nearby  C u  moments,  thus  confirming  static 
magnetic  order.  In  contrast,  the  time  spectra  of  supercon¬ 
ducting  specimens,  Y Ba^Cu^Oe.s*  and  Aai.gsFrg  ,7iCuO^, 
showed  only  very  slow  depolarization,  wholse  rate  can  be 
accounted  for  by  nuclear  dipolar  fields.  This  rules  out  the 
existence  of  static  magnetic  order  (with  the  accuracy  better 
than  0.01  ms  per  Cu)  in  good  superconducting  specimens. 
This  demonstrates  that  static  magnetic  order  is  not  an  in¬ 
trinsic  property  of  superconducting  hole-doped  high-Tc  sys¬ 
tems.  The  results  of  ref.  7  may  be  due  to  inhomogeneity  of 
the  specimen. 

Magnetic  order  —  electron-doped  systems 

We  have  also  performed9  mSR  studies  on  various  spec¬ 
imens  of  the  electron-doped  system  Nd2-zCezCuOt  (0.0  < 
i  <  0.18;  reduced  and  non-reduced).  Figure  2  shows  typical 
zero-field  pSR  spectra.  The  pure  parent  system  NdjCuOi 
exhibits  static  magnetic  order  with  Tpt  =  250A’.  We  have 
also  found  a  signature  of  spin  re-orientation  at  T  ~  80A'  and 
T  ~  30 A'  in  this  system,  due  to  the  instability  of  spin  direc¬ 
tions  in  the  tetragonal  crystal  structure.  The  fast  spin  deoo- 
larization  in  non-superconducting  N  d\.gCe0,iCuOt  confirms 
that  the  static  magnetic  order  extends  in  the  doped  system. 
In  contrast,  the  superconducting  specimen  with  x  =  0.16, 
reduced,  does  not  show  any  signature  of  static  order  of  Cu 
spins  above  T  —  5 A'. 
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Fig.  2.  Zero-field  p SR  spectra  in  the  electron-doped  system 
Nd7-tCeTCuO<  with  x  =  0.0  and 0.1  (non- superconducting), 
and  x  -  0.1G  (superconducting)  reduced  ceramic  specimens 
(from  Luke  el  al.9). 
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Fig.  3.  Magnetic  and  superconducting  phase  diagram  of 
Ndi_zCezCuOi .  Magnetic  phase  boundary  us  determined 
by  pSR,  while  superconducting  phase  boundary  by  resistivity 
(from  Luke  el  al.9). 


Based  on  these  results,  we  obtained  the  magnetic  phase 
diagram  as  shown  in  Fig.  3.  Similarly  to  hole-doped  sytems. 
the  static  magnetic  order  disappears  around  the  concentra¬ 
tion  where  the  superconductivity  appears.  Here  again,  there 
is  no  clear  evidence  for  the  intrinsic  coexistence  of  supercon¬ 
ductivity  and  the  magnetic  order  of  Cu  spins.  In  the  electron- 
doped  systems,  the  magnetic  order  extends  to  much  higher 
impurity  concentration  (i  ~  0.15)  compared  to  the  case  of 
the  hole-doped  214  system.  A  hole  doped  in  the  CuOi  nlane 
is  believed  to  stay  at  the  oxygen  site  as  Ox~ ,  creating  a  frus¬ 
trated  bond  within  the  Cu  T  07~Cu  1  antiferromagnetic  net 
work.  The  electron  carrier,  in  contrast,  goes  to  the  Cu  atom, 
making  it  as  non  magnetic  Cu  Tl-  Thus  the  electron  doping 
is  equivalent  to  dilution  of  the  magnetic  moments,  which  has 
much  milder  effect  in  destroying  spin  order  compared  to  frus¬ 
tration.  This  argument  qualitatively  explains  the  extended 
region  of  the  antiferromagnetic  phase  in  Fig.  3. 

Penetration  depth  -  Tc  v.s.  n ,/m' 

To  measure  the  magnetic  field  penetration  depth  A,  we 
apply  transverse  external  magnetic  field  Hett  (Hci  <  Hez,  < 
Hcl\  typically  HeI,  =  2  -  20 kG),  and  observe  the  depolariza¬ 
tion  of  the  muon  spin  precession:  for  simplicity10,  one  usu¬ 
ally  assumes  the  Gaussian  relaxation  function  exp(— o2f  (2) 
to  define  the  relaxation  rate  o.  Figure  4  shows  examples  of 
the  tempeiature  dependence  o(T)  observed  in  two  different 
high-Tc  superconductors,  measured  by  cooling  in  Htzt  (field 
cooling).  In  the  superconducting  state,  Hezt  penetrates  the 
specimen,  forming  a  lattice  of  fiux  vortices.  The  local  field 
B  becomes  inhomogeneous:  B  is  largest  at  the  core,  while 
smallest  in  the  midpoint  beween  adjacent  fiux  cores.  As  an 
ensamble  average  signal  from  many  muons  stopped  at  differ¬ 
ent  locations  with  respect  to  the  fiux  lattice,  pSR  samples  the 
field  distribution  in  the  vortex  state.  The  London  equation 
leads  to  AB  <x  1/A2,  where  A5  =  (<  ( B -  <  B  >)'  >)1/'- 


Therefore,  by  measuring  the  relaxation  rate  a ,  one  can  de¬ 
termine  the  penetration  depth,  which  is  related  to  the  super¬ 
conducting  carrier  density  n,  divided  by  the  effective  mass 
m'  as 

a  a  AB  oc  1/A5  <x  n,/m‘. 

We  performed11  the  measurements  of  o{T )  on  more  than 
20  specimens  (unoriented,  sintered  disks)  of  various  high-Tc 
superconductors.  Assuming  the  angular  averaging  process  for 
anisotropy  of  A  to  be  the  same  for  different  specimens,  one 
can  regard  a(T  — >  0)  to  represent  the  ground  state  values  of 
njnt'  of  each  system.  We  selected  the  results  from  speci¬ 
mens  which  show  sharp  changes  of  a  at  Tc  with  greater  than 
a  90  %  superconducting  volume  fraction.  Figure  5  shows  Tc 
and  o(T  — *  0),  both  determined  by  fiSR.  With  increasing  hole 
doping,  Te  initially  increases,  following  a  linear  relation  with 
n,/m*  universal  to  systems  with  single,  double,  and  triple 
CuOj  layers.  Then,  Tc  shows  saturation  and  suopression  in 
the  heavily-doped  region.  Recently,  we  have  added15  two 
points  on  (Y\^IPrI)BaiCu30r  with  x=0.15  and  0.3  in  Fig. 
4,  which  agree  quite  well  with  the  results  on  Y BaiCu3Or. 
This  suggests  that  Tc  is  determined  solely  by  the  carrier  den¬ 
sity,  regardless  of  the  source  of  the  doping. 

;iSR  is  unique  in  providing  information  directly  related 
to  the  superconducting  carrier  density.  The  results  in  Fig. 
5  should  be  explained  by  successful  theory  of  high-Tc  super¬ 
conductivity.  The  horizontal  axis  o  <x  n,/m‘  is  proportional 
to  Fermi  energy  ej  of  non-interacting  2-dimensional  electron 
gas.  Therefore,  one  possible  interpretation  is  to  view  the 
straight  line  as  indicating  a  </,  which  is  not  expected  in 
the  weak-coupling  BCS  model  but  is  conceivable  when  the 
energy  scale  of  the  pairing  interaction  is  larger  than  ej  (i.e., 
no  retardation13 ). 
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Fig.  4 ■  The  temperature  dependence  if  the  muon  spin  relax¬ 
ation  rate  <7  observed  in  superconducting  Tl-Pb-Sr-Ca-Cu-0 
systems  (with  double  and  triple  Cu03  layers J  in  a  transverse 
external  magnetic  field  of  2.5  kG  (from  Uemura  et  al. 
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Fig.  5.  The  superconducting  transition  temperature  Tc  plot¬ 
ted  versus  low-temeprature  muon  spin  relaxation  rate  o(T  — 
0)  of  various  planar  hole-doped  high-Tc  cuprates  (see  rtf.  11 
for  details).  The  closed  square  points  represent 
(Yi-tPrJB^CuiOr  With  x  =  0.15  and  0.3.  (from  Uemura 
el  al.11). 
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Fig.  6.  The  temperature  dependence  of  the  magnetic 

field  penetration  depth  Aj|  derived  from  the  muon  spin  relax¬ 
ation  rate  o  in  a  c-axis  oriented  sintered  ceramic  specimen 
ofYBaiCusOi  with  the  external  field  HeIi  ~  13fcC  applied 
imrallel  to  the  c-axts.  The  solid  line  represents  the  fit  to  the 
two-fluid  model,  (from  Uemura  et  al)  )■ 


Penetration  depth  —  energy  gap 

Thermai  excitations  across  the  energy  gap  destroy  the 
Cooper  pairs,  reduce  n,,  and  make  the  penetration  depth  at 
finite  temperatures  longer  than  A (T  =  0).  Figure  6  shows 
A (T)  Bdetermined  in  a  c-axis  oriented  ceramic  specimen  of 
YBa7Cu307  with  Hczt  ||  c-axis1.  A(T)  agrees  reasonably 
well  with  the  two  fluid  model,  indicating  that  the  energy  gap 
is  finite  without  any  detectable  anomalous  zeros.  Essentially 
the  same  results  are  obtained  on  non-oriented  specimens  of 
all  the  other  hole-doped  planer  systems  (214,  2212,  123,  2223, 
and  similar  systems).  Experiments  are  now  underway  to  ex¬ 
tend  the  pSR  measurements  of  A(T)  to  heavy-fermion  and 
organic  superconductors. 

Flux  pinning 

Recently,  we  have  found1*  that  pSR  can  also  detect  flux 
pinning  phenomena.  See  Fig.  7.  When  an  external  field  is 
applioed  above  Tc  followed  by  field  cooling  (FC),  there  is  no 
anomaly  in  the  er(T)  curve.  With  Hett  applied  perpendic¬ 
ular  to  the  face  of  the  thin  disk  specimen,  the  total  mag¬ 
netic  induction  B  is  the  same  above  and  below  Tc.  Then, 
the  flux  vortex  does  not  have  to  move  any  macroscopic  dis¬ 
tance.  and  the  FC  pSR  results  are  insensitive  to  flux  pinning. 
In  contrast,  if  one  applies  the  external  field  after  cooling  in 
zero  field  (ZFC),  the  flux  vortices  must  move  macroscopic 
distances.  If  the  flux  is  trapped  during  this  motion,  the  re¬ 
sulting  field  distribution  would  be  more  inhomogeneous  than 
the  case  of  FC.  Therefore,  the  difference  between  the  FC  and 
ZFC  data  in  Fig.  7  shows  that  the  flux  pinning  takes  place 
below  T  ~  30 A'  in  the  Bi  2212  system.  This  pinning  tem¬ 
perature  agrees  very  well  with  those  found  in  resistivity  and 
mechanical  measurements15.  This  serves  to  confirm  that  the 
phenomena  observed  in  these  bulk  measurements  were  indeed 
related  to  the  motion  of  the  flux  vortex. 

Sintered  Bi*Sr,CaCujO» 
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Fig,  7.  Comparison  of  the  muon  spin  relaxation  rate  o  in  a 
sintered  ceramic  sj>ecimen  of  Bi3SriCaCu30t  measured  with 
the  field  cooling  (FC)  and  zero-field  cooling  (ZFC)  procedures 
with  the  external  field  of  f.5  kC.  The  difference  between  FC 
and  ZFC  results  below  T  =  30  h'  is  due  to  flux  pinning,  (from 
Sternheb  el  at. M ). 


Conclusions 

The  temperature/concentration  phase  diagrams  of  pla¬ 
nar  cuprate  high-7"c  systems,  obtained  by  pSK  both  in  hole- 
doped  and  electron-doped  compounds,  indicate  the  existence 
of  critical  carrier  concentrations  which  separate  the  semi¬ 
conducting  region  with  statically  ordered  Cu  moments  and 
the  superconducting  region  with  metallic  normal-state  con¬ 
ductivity.  The  boundary  between  these  two  regions  is  very- 
sharp,  and  there  is  no  signature  of  intrinsic  coexistence.  This 
strongly  suggests  that  the  mechanism  of  high-Tc  supercon¬ 
ductivity  is  intimately  related  to  magnetism.  The  univer¬ 
sal  correlation  between  n,/m *  and  Tc  in  the  superconduct¬ 
ing  state,  as  well  as  the  existance  of  finite  energy  gap  with¬ 
out  anomalous  zeros,  will  be  useful  for  the  determination  of 
pairing  and  condensation  mechanisms.  Ongoing  pSR  exper¬ 
iments  will  continue  to  provide  information  on  these  aspects 
as  well  as  on  flux  pinning  phenomena  in  various  supercon¬ 
ducting  systems. 
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The  crystal  chemistry  copper  oxide  superconductors  is  reviewed  with  respect  to  defects,  stability  and 
the  n-type  phases.  It  is  concluded  that  the  n-type  copper  oxide  superconductors  are  genuine  exam¬ 
ples  of  reduced  systems  which  any  viable  theory  for  high  Tc  must  explain.  Defects  are  crucial  for  the 
control  of  carrier  concentration  and  presumably  critical  current.  They  also  cause  troublesome  and 
pervasive  inhomogeneities.  Most,  and  possibly  all,  high  Tc  superconductors  are  thermodynamically 
unstable  at  room  temperature  and  below. 


1.  INTRODUCTION 

I  have  recently  reviewed  the  crystal  chemistry 
of  oxide  superconductors1'3;  thus,  this  paper  will 
not  be  yet  another  overview  of  the  crystal  chem¬ 
istry  of  these  materials.  Instead,  I  focus  on  the 
crystal  chemistry  issues  relating  to  defects,  stabil¬ 
ity  and  the  n-type  copper  oxides.  The  role  of 
defects  is  critical  because  they  control  carrier  con¬ 
centration,  inhomogeneities  and  presumably  flux 
pinning.  The  apparent  lack  of  thermodynamic 
stability  for  high  Tc  superconductors  is  a  crucial 
factor  in  the  synthesis  of  these  materials  and 
must  also  be  considered  for  any  theory  of  the 
mechanism  for  high  Tc. 

2.  n-TYPE  COPPER  OXIDES 

From  the  point  of  view  of  our  understanding  of 
copper  oxide  based  superconductors,  the  n-type 
copper  oxides  found  by  Tokura  et  al*  represent  the 
most  significant  new  superconductors  since  the 
discovery  of  Bednorz  and  Muller.  Different  mech¬ 
anisms  for  n-  and  p-type  materials  is  a  remote 
possibility.  Thus,  an  important  part  of  the  correct 
theory  must  be  that  p-  and  n-type  doping  give 
essentially  the  same  effect.  This  then  eliminates 
from  consideration  those  theories  dependent  on 
some  special  feature  of  an  "oxygen  2p  band",  espe¬ 
cially  a  re*  oxygen  2p  band. 

Evidence  for  the  n-type  character  of  certain 
copper  oxides  does  not  come  merely  from  electri¬ 
cal  transport  data.  There  is  also  overwhelming 
evidence  from  the  chemistry  side.  Several  papers 
presented  at  this  conference5  8  have  amply 
demonstrated  that  these  materials  are  now  as 
well  characterized  as  are  most  of  the  p-type  copper 
oxides. 

The  materials  first  discovered  to  be  n-type  cop¬ 
per  oxides  are  the  R2  xCexCu04±y  series  where  R 
may  be  Pr,  Nd,  Srn,  Eu  or  Gd.*  Initially,  some 
doubts  were  expressed  about  the  description  of 
these  phases  because  cerium  is  well  known  in 


both  the  tetravalent  and  trivalent  oxidation  states. 
This  issue  was  clouded  by  widespread  confusion 
concerning  the  meaning  of  oxidation  states,  ex¬ 
emplified  by  the  absurd  notion  the  cerium  might 
not  be  tetravalent  in  CeC>2.  Spectroscopic  studies 
have  now  shown  that  cerium  is  tetravalent  in  the 
R2-xCexCu04  superconductors.  Actually,  the 
oxidation  state  of  Ce  was  already  clear  from 
examination  of  the  unit  cell  dimensions  of  the 
R2-xCexCu04  phases.  We  know  that  Ce 17  is 
smaller  than  the  R111  cations  and  that  Ce111  is 
larger  than  the  R111  cations  being  used.  Thus,  we 
expect  to  see  the  c  axis  increase  if  cerium  were 
trivalent.  In  fact,  the  c  axis  decreases  with 
cerium  substitution  showing  that  the  cerium  is 
tetravalent.  The  a  axis  slightly  increases  with  Ce 
substitution  despite  the  smaller  size  of  Ceiv.  This 
occurs  because  the  a  axis  dimension  is  dominated 
by  the  Cu-0  distance  of  the  highly  covalent  Cu-0 
bond.  The  highest  filled  levels  in  a  Cu11  oxide  are 
strongly  antibonding.  Thus,  removing  electrons 
causes  a  decrease  in  the  Cu-0  bond  length. 

Adding  electrons  results  in  an  increase  of  the 
bond  distance  and  thus  an  increase  in  the  a  cell 
dimension. 

In  the  R2.*Th*Cu04  and  Nd2Cu04.»F*  phases,5 
there  are  no  oxidation  state  issues  for  the  dopants. 
The  remaining  issue  is  the  question  of  precise 
oxygen  content.  Tarascon  presented  very  careful 
studies  at  this  conference5  showing  that  super¬ 
conducting  properties  were  optimized  when  both 
oxygen  interstitials  and  vacancies  were  mini¬ 
mized.  Thus  from  a  chemistry  point  of  view,  these 
materials  are  dearly  reduced,  n-type  systems. 

This  is  confirmed  by  chemical  titrations  for  oxida¬ 
tion  states.  Various  spectroscopic  studies  have 
also  confirmed  the  reduced  nature  of  copper  in 
these  systems. 

Hall  and  Seebeck  measurements  have  shown 
dominantly  n-type  behavior  for  these  materials. 


Tin*  hub  been  true  lor  both  ceramic  and  single 
crystal  materials,  regardless  of  whether  the 
measurement  is  taken  parallel  or  perpendicular 
c  axis.  There  are  conditions  for  which  positive 
Hall  coefficients  have  been  found  for  the 
R2 -*Ce„Cu04  phases,  but  there  are  also  cases 
where  negative  Hall  coefficients  are  found  for 
YBa2Cu3Ofl.n1  phases. 

The  coordination  of  Cu1  to  oxygen  is  normally 
two-fold  linear.  Thus,  Cu1  in  the  square  planar 
coordination  of  the  T  structure  might  seem  im¬ 
plausible.  However,  there  are  examples  where 
the  coordination  of  Cu1  to  oxygen  is  four  fold  and 
not  tetrahedral.9  More  importantly,  examination 
of  the  causes  of  two-fold  linear  coordination  for 
Cu1  leads  one  to  believe  that  square  planar  coordi¬ 
nation  for  Cu1  is  more  favorable  in  oxides  than  is 
tetrahedral  coordination.  Frequently,  the  bonding 
of  Cu1  is  simply  viewed  as  utilizing  s  and  p 
orbitals  of  copper  because  the  3d  shell  is  filled. 
However,  the  pervasiveness  of  "sp"  bonding  and 
lack  of  sp2  and  sp3  bonding  for  Cu1  in  oxides  can¬ 
not  be  understood  on  this  basis.  We  must  assume 
that  the  filled  3d  core  interferes  with  s-p  bonding. 
Thus,  the  preiequisite  for  bond  formation  is  polar¬ 
ization  of  the  filled  core  through  s-d  hybridiza¬ 
tion.10  The  usual  occurrence  is  that  the  core 
polarizes  from  a  sphere  to  a  disk.  Thus,  only  two 
strong  bonds  may  form  along  the  unique  axis  of 
this  disk.  However,  s-d  hybridization  could  as 
well  polarize  the  core  to  give  an  egg-shaped  ellip¬ 
soid.  Then,  the  strong  Cu^O  bonds  would  form 
perpendicular  to  the  unique  axis,  and  this  is 
consistent  with  square  planar  coordination  for 
Cu1.  This  situation  is  analogous  to  the  Jahn 
Teller  effect  found  for  d8  cations  such  as  Cu111. 

The  usual  occurrence  of  this  Jahn  Teller  distor¬ 
tion  of  an  octahedron  is  that  we  see  four  short  and 
two  long  bonds.  A  two-short,  four-long  configura¬ 
tion  is  energetically  equivalent  to  the  first  approx¬ 
imation.  However,  this  configuration  is  only 
rarely  found.  Likewise,  four  short  bonds  for  Cu1  is 
only  rarely  found  in  oxides  even  though  its  ener¬ 
getics  are  not  so  unfavorable  in  the  square  planar 
configuration. 

The  existence  of  both  n-  and  p-type  copper 
oxides  suggests  a  schematic  phase  diagram  and 
band  structure  as  shown  in  Figure  1.  The  sym¬ 
metry  is  exaggerated  in  both  parts  of  this  figure, 
but  I  doubt  that  the  deviations  from  this  ideal 
symmetry  are  important  for  a  basic  understand¬ 
ing  of  these  materials.  At  the  midpoint  of  the 
schematic  phase  diagram,  all  copper  is  divalent 
and  the  o*  band  is  split  to  give  an  antiferromag¬ 
netic  insulator.  States  near  the  Fermi  level  are 
derived  nearly  equally  from  oxygen  and  copper; 
thus,  labeling  bands  near  the  Fermi  level  as  Cu  3d 
or  as  O  2p  is  misleading. 


FIGURE  1 

Schematic  phase  diagram  (lower)  and  schematic 
band  structures  (upper)  for  A2CUO4  phases,  both 
n-  and  p-types 

Reduction  of  Cu11,  adding  electrons  from  the 
midpoint,  leads  to  n-type  materials.  Oxidation  of 
Cu1*,  removing  of  electrons  from  midpoint,  leads 
to  p-type  materials.  Doping  either  way  destroys 
the  antiferromagnetic  state.  Sufficient  doping 
produces  conductors  which  superconduct.  The 
materials  with  the  highest  doping  levels  are  in¬ 
completely  characterized.  They  are  difficult  to 
prepare  single  phase,  and  single  crystals  are  not 
yet  available.  Nonetheless,  current  evidence  sug- 
ests  that  the  carrier  type  switches  as  supercon- 
uctivity  disappears.  This  is  consistent  with  a 
view  that  the  splitting  of  the  o*  band  disappears  at 
high  doping  levels. 

The  x*  band,  which  overlaps  the  filled  part  of 
the  a*  band!,  is  ignored  in  this  figure.  The  observ¬ 
ed  symmetry  between  electron  and  hole  carriers 
suggests  that  this  x*  band  can  be  safely  ignored 
for  a  basic  understanding  of  these  materials. 
However,  until  Te 's  for  n-type  materials  match 
those  for  p-type  materials,  it  is  natural  to  think 
that  the  x*  band  plays  some  role  in  enhancing  Tc. 

The  symmetry  of  the  phase  diagram  (Fig.  1) 
has  been  exaggerated.  Nonetheless,  the  optimum 
doping  level  is  essentially  the  same  for  the  p-  and 
n-doping,  i.e.  about  0.15  electrons  or  holes  per  cop¬ 
per.  This  is  additional  evidence  that  reference  to 
copper  oxidation  states  is  appropriate  despite 


the  high  covalcncy  with  strong  copper-oxygen 
admixture  in  the  region  of  the  Fermi  level. 

There  are  several  significant  departures  from 
the  symmetry  shown  in  the  schematic  phase  dia¬ 
gram  (Fig.  1).  For  the  n-type  materials,  the  mag¬ 
netic  insulating  state  actually  disappears  more 
slowly. 1 1  The  superconducting  state  is  observed 
over  a  more  narrow  range  of  doping,  and  Tc 
reaches  only  about  25  K  instead  of  about  35  K  for 
the  p-type  phases.  This  lack  of  symmetry  may 
well  be  related  to  the  fact  that  the  o*  states  gener¬ 
ally  derive  more  from  copper  and  less  from  oxygen 
as  the  cr*  band  fills. 

The  scale  has  been  purposely  left  off  the  hori¬ 
zontal  axis  of  the  phase  diagram.  This  is  because 
all  doped  materials  are  intrinsically  inhomoge¬ 
neous,  and  the  scale  of  the  horizontal  axis  thus 
becomes  confused.  We  will  come  back  to  this  in¬ 
homogeneity  issue  in  the  next  section.  However, 
this  confusion  in  no  way  invalidates  the  general 
trends  shown  in  the  phase  diagram  section  of 
Figure  1. 

One  of  the  most  inappropriate  questions  ever 
asked  of  the  p-type  copper  oxides  was  whether  the 
hole  was  an  oxygen  or  copper.  The  question  was 
inappropriate  because  of  the  high  degree  of  cop¬ 
per-oxygen  admixture  near  the  Fermi  level.  Also, 
we  know  that  the  holes  would  not  develop  reason¬ 
able  mobility  if  they  did  not  spend  significant  time 
on  both  copper  and  oxygen.  Actually,  it  is  best  to 
think  of  the  holes  as  being  on  neither  copper  nor 
oxygen.  Instead,  they  reside  in  the  bonds  between 
copper  and  oxygen. 

Once  the  question  was  asked  as  to  whether  the 
holes  were  on  copper  or  oxygen,  people  rose  to  the 
challenge  to  provide  an  answer.  Evidence  was 
found  for  oxygen  2p  character  for  the  holes.  Thus, 
some  too  quickly  concluded  that  the  holes  were  on 
oxygen  alone.  However,  equally  good  evidence  has 
been  found  that  the  holes  have  copper  3d  charac¬ 
ter.  This  evidence  come  from  XAS  studies12  as 
well  as  XPS  studies  of  both  core  peaks13  and  the 
valence  region.14  The  best  current  estimate  of  the 
relative  contributions  of  oxygen  and  copper  at  the 
Fermi  level  is  that  the  states  derived  about  70% 
from  oxygen  and  30%  from  copper. 14  Taking  into 
account  the  composition  of  the  sheets,  CuC>2,  the 
conclusion  is  that  the  individual  atoms  of  the 
sheets  are  contributing  nearly  equally  to  states 
near  the  Fermi  level. 

3.  DEFECTS 

Defects  in  the  copper  oxide  systems  are  both  a 
curse  and  a  blessing.  We  could  not  control  carrier 
concentrations  without  the  defects,  and  it  is  varia¬ 
tion  of  carrier  concentration  that  has  given  us  our 
most  interesting  trends  with  Te,  Tn,  etc.  Defects 
are  also  the  hope  for  improving  the  critical  cur¬ 


rent  in  the  copper  oxido  based  superconductors. 
Finally,  defects  cause  inhomogonoities  which 
plague  us  continuously. 

Solid  solutions,  including  systems  of  the  type 
Laa.xSr.CuO*.  are  intrinsically  inhomogeneous. 
The  substitution  of  Sr  for  La  creates  a  point  defect. 
Electronic  properties  in  the  vicinity  of  Sr  are 
different  than  elsewhere  in  the  La2Cu04  matrix. 

If  one  uses  the  coherence  length  to  set  the  scale  for 
the  size  of  significant  inhomogeneity,  we  find  that 
there  is  no  value  of  x  for  which  the  inhomo- 
geneities  can  be  safely  ignored.  Flux  exclusion 
data  indicate  that  the  most  homogeneous  super¬ 
conductor  in  the  La2-xSrxCu04  system  is  for  x  of 
about  0.15.  However,  even  this  composition  cannot 
be  regarded  as  homogeneous  because  a  random 
distribution  of  Sr  at  this  concentration  in  the 
La2Cu04  matrix  would  give  significant  Sr-free 
regions  which  are  larger  than  twice  the  coherence 
length. 

Inhomogeneities  should  not  be  mistaken  for 
indications  of  phase  separation.  A  material  may 
be  strictly  single  phase  in  the  usual  thermody¬ 
namic  sense  and  yet  be  very  inhomogeneous  in  an 
electronic  sense.  The  basic  problem  is  that  solid 
solutions  and  doped  systems  are  intrinsically 
inhomogeneous.  Sometimes,  the  inhomogenei¬ 
ties  are  not  important  for  the  properties  of  in¬ 
terest.  However,  for  the  oxide  superconductors 
such  as  those  in  the  La2-xSrxCuC>4  system,  these 
inhomogeneities  cannot  be  ignored  with  respect  to 
either  magnetic  or  electronic  properties.  One  of 
the  disturbing  consequences  of  this  fact  is  that 
some  of  our  probes  will  invariably  be  influenced 
more  by  the  low  x  regions  and  other  probes  will  be 
more  influenced  by  the  high  x  region  of  our  single 
phase  solid  solutions.  Thus,  those  performing 
different  measurements  may  have  great  difficulty 
reconciling  their  results  because  they  are  probing 
different  regions  of  a  heterogeneous,  single-phase 
material. 

There  are  also  several  extended  defects  to 
consider  in  the  oxide  superconductors.  In  the 
YBa2Cu307  type  phases,  there  are  both  1 10  twin¬ 
ning  and  90°  misorientation  defects.  The  latter  is 
also  sometimes  referred  to  as  twinning.  However, 
the  110  twinning  and  the  90°  misorientations  are 
distinctly  different.  The  former,  coherent  twin¬ 
ning,  relieves  lattice  strain  while  the  latter,  inco¬ 
herent  twinning,  causes  lattice  strain.  There  has 
been  much  speculation  about  the  roles  of  these  two 
extended  defects  in  YBa2Cu3<I>7  phases.  Another 
common  extended  defect  is  stacking  faults  or 
intergrowth.  These  are  particularly  common  in 
the  Tl,Bi/Ba,Sr/Ca,Cu/0  superconductors. 

A  new  extended  defect  has  recently  been  dis¬ 
covered  by  Gai15  in  the  L42CUO4  structure.  This 
defect  is  essentially  absent  in  La2CuC>4  itself. 


However  in  the  La2-iSr,,Cu04  system,  this  defect 
increases  with  x.  at  least  up  to  x  =  0. 15.  The  defect 
may  be  described  as  a  shear,  and  oxygen  is  elimi¬ 
nated  at  this  defect.  Thus,  this  system  is  more 
properly  regarded  as  La2.xSrxCu04.y  where  the 
extended  defect  and  y  tend  to  increase  with  in¬ 
creasing  x.  Thus,  the  carrier  concentration  in 
these  phases  is  expected  to  be  somewhat  less  than 
that  predicted  by  the  value  of  x  alone.  However, 
this  effect  could  be  balanced  off  to  some  extent  by 
interstitial  oxygen. 

It  has  been  mysterious  why  Tc  drops  when 
rare  earths  such  as  Nd  substituted  for  La  in 
superconductors  of  the  type  La2.xSrxCu04. 
Appreciable  composition  ranges  of  the  type  La2-x- 
ySr,RyCu04  may  be  prepared  where  R  is  Pr,  Nd, 
Sm,  Eu  or  Gd.16  The  Cu-0  distance  in  the  Cu02 
sheets  decreases  as  R  becomes  smaller  and  as  y 
becomes  larger.  Nonetheless,  the  tendency  is  for 
Tc  to  drop  with  this  substitution,  the  opposite  of  the 
usual  correlation  between  Tc  and  Cu-0  distance. 
We  decided  that  a  missing  piece  of  this  puzzle  was 
the  chemical  analysis  for  the  oxidation  state  of 
copper.  On  performing  this  analysis  on  many 
samples  of  the  type  La2-x-ySrxRyCu04,  we  found 
that  for  fixed  x  the  Cu111  content  tends  to  decrease 
as  y  increases  and  as  R  becomes  smaller.17  In 
fact,  these  materials  have  the  very  familiar  trend 
of  Tc  with  Cu111  concentration.  This  strong  correl¬ 
ation  obscures  the  possible  correlation  of  TV  with 
Cu-0  distance.  The  shear  type  defects  found  by 
Gai  in  La2-xSrxCu04  are  also  pronounced  in 
La2x-ySrxRyCu04  phases.17  We  may  speculate 
that  these  defects  are  now  more  numerous  and 
that  this  has  resulted  in  a  greater  oxygen 
deficiency. 

Not  all  oxide  superconductors  depend  on 
defects  for  all  the  important  Cu111  or  Cu1  content. 
The  ideal  formulae  for  YBa2Cu3C>7,  YJ^CusOs 
and  Y2BaiCu70is  are  all  apparently  appropriate 
to  give  high  Tc  superconductors.  These  then  are 
the  materials  of  choice  if  one  wishes  to  study 
homogeneous  materials.  However,  these  materi¬ 
als  are  prone  to  have  defects  such  as  twins, 
misorientations  and  intergrowths.  Therefore, 
their  defect  concentration  may  be  significant 
Another  problem  with  all  of  these  phases  is  that 
much  of  the  copper  is  not  in  the  Cu(>2  sheets.  This 
is  not  true  of  the  other  copper  oxide  based  super¬ 
conductors  where  the  carrier  concentration  is 
controlled  through  the  point  defect  concentration. 

Band  structure  calculations  had  suggested 
that  the  carrier  concentration  in  BfeSrgCan-i 
Cun04*2n  phases  is  controlled  by  overlap  of  the  Bi 
6s  band  with  the  Fermi  level.  However,  this  con¬ 
clusion  is  very  much  at  odds  with  reasonable 
chemistry1  and  apparently  resulted  from  the 


highly  idealized  structure  used  for  tho  calcula¬ 
tions.  The  carrier  concentration  in  these  phases 
is  instead  due  to  deviations  from  the  ideal  stoi¬ 
chiometry,  i.e.  point  defects.2 
4.  STABILITY 

There  is  good  reason  to  believe  that  most  (or  all) 
of  the  high  Tc  superconductors  are  not  thermody¬ 
namically  stable  at  room  temperature  and  below. 
This  is  not  necessarily  a  concern  in  practical 
terms.  After  all,  diamond  and  window  glass  are 
also  mctastable  under  ambient  conditions. 
However,  optimum  synthesis  approaches  to 
mctastable  materials  are  very  different  than  those 
employed  for  thermodynamically  stable  phases. 
Also,  it  would  seem  that  a  viable  theory  must 
suggest  an  unstable  lattice  in  a  high  temperature 
superconductor.  This  was,  of  course,  relatively 
straightforward  using  classical  BCS  theory. 

The  standard  synthesis  of  YBa2Cu306+x 
phases  starts  under  conditions  where  x  in  the 
product  is  so  small  (x  -  0.3)  that  superconducting 
materials  are  not  obtained.  These  phases  must  be 
oxidized  to  obtain  a  high  Tc  superconductor.  The 
oxidation  of  YBa2Cu3C>6  to  YBa2Cua07  is  an  inter¬ 
calation  reaction.  Products  made  by  intercalation 
reactions  are  nearly  always  metastable.  Thus, 
they  cannot  generally  be  made  in  a  direct,  one-step 
process.  All  attempts  to  prepare  YBa2Cu3C>7  in  a 
one  step  process  have  failed.  Serious  attempts  to 
obtain  true  equilibrium  for  YBa2Cu3C>6+x  phases 
indicated  that  the  highest  value  of  x  for  a  thermo¬ 
dynamically  stable  phase  is  about  0.7. 18  However, 
even  this  phase  is  most  likely  not  stable  at  room 
temperature  and  below.  After  all,  the  most  stable 
phase  in  the  YBa2Cu3O6.n1  series  (x  =  0.0)  is 
apparently  not  stable  below  about  500°C. 

The  structures  of  YBa2Cu40g  and  Y2Ba4Cu70is 
are,  of  course,  very  closely  related  to  that  of 
YBa2Cu307,  and  it  would  appear  that  both 
YBa2Cu408  and  Y2Bo4Cu70xs  are  thermodynami¬ 
cally  stable  phases  under  the  conditions  where 
they  are  prepared.  There  is  no  current  evidence 
that  these  materials  become  unstable  at  lower 
temperature,  very  unlike  the  next  compounds  to 
be  considered. 

Our  highest  Tc's  have  been  found  for  com¬ 
pounds  which  are  ideally  represented  as 
(A'0)xA2Can.iCun02n+2  where  A'  is  T1  or  Bi  and 
A  is  Ba  or  Sr.  Actually,  these  phases  are  always 
ridden  with  defects  resulting  in  substantial  devia¬ 
tions  from  this  ideal  formulation.  These  phases 
have  only  been  prepared  at  high  temperatures, 
very  close  to  melting  temperatures.  All  low 
temperature  synthesis  routes  have  failed.  Also, 
exotherms  of  compound  formation  have  never 
been  noted  for  these  compounds  when  they  are 
prepared  from  the  binary  oxides.  Thus,  there  is 
good  reason  to  believe  that  thrxo  phases  ore 


REFERENCES 


endothermic  compounds,  i.e.,  they  are  stabilizod 
by  entropy  and  are  not  stable  at  lower  tempera¬ 
tures. 

Entropy  stabilized  oxide  compounds  are  rare 
but  known.  Well  studied  examples  are  CuFe204. 
Fe^TiOa  and  AleSi20i3.19  In  all  of  these  cases,  the 
disorder  giving  the  entropy  for  compound  stabi¬ 
lization  is  in  the  form  of  point  defects. 

Thus,  the  (A'0)*A2Can.iCun02n+2  supercon¬ 
ductors  appear  to  cases  of  entropy  stabilized 
phases,  i.e.,  they  are  stable  under  their  synthesis 
conditions  but  not  at  lower  temperatures.  They 
are  defect  ridden,  and  apparently  these  defects  are 
required  for  the  stabilization  of  these  phases.  This 
also  means  that  these  phases  will  tend  to  be 
inhomogeneous  even  though  we  frequently  repre¬ 
sent  them  with  ideal  formulae. 

Solid  solutions  or  doped  systems  are  trivial 
examples  of  entropy  stabilized  systems.  Thus, 
superconductors  such  as  La2-xSrxCuO<  and 
Nd2-*Ce*Cu04  are  not  expected  to  be  thermody¬ 
namically  stable  at  low  temperatures.  In  fact,  the 
thermodynamic  stability  range  appears  to  be  suf¬ 
ficiently  restricted  to  preclude  low  temperature 
synthesis  of  these  phases. 

5.  CONCLUSIONS 

The  chemistry  case  for  the  n-type  copper  oxides 
is  overwhelming.  These  materials  are  what  they 
appear  to  be,  and  successful  theory  must  account 
for  this  symmetry  at  the  Fermi  level. 

Defects  control  many  of  the  properties  of  inter¬ 
est.  We  are  still  struggling  to  control  and  mea¬ 
sure  defect  concentration.  Much  work  remains  to 
relate  defects  to  flux  pinning.  The  defects  lead 
inescapably  to  inhomogeneities  in  most  of  the 
oxide  superconductors.  We  cannot  afford  to 
ignore  these  inhomogeneous  materials  which 
make  up  the  bulk  of  the  high  temperature  super¬ 
conductors.  Instead,  we  must  learn  to  model  this 
inhomogeneity. 

Most  oxide  superconductors  are  thermody¬ 
namically  unstable.  It  remains  to  be  shown 
whether  or  not  they  are  all  metastable.  For  those 
we  now  know  to  be  metastable,  rational  synthesis 
techniques  must  take  into  account  this  lack  of 
thermodynamic  stability. 
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